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Abstract 

This  thesis  describes  the  application  of  the  multiple-input  multiple-ouput  (MIMO)  Quantita¬ 
tive  Feedback  Theory  (QFT)  design  technique  to  the  design  of  a  digital  flight  control  system  for  the 
Lambda  Unmanned  Research  Vehicle  (URV).  The  QFT  technique  allows  the  synthesis  of  a  control 
system  which  is  robust  in  the  presence  of  structured  plant  uncertainties.  Uncertainties  considered 
in  this  design  are  the  aircraft’s  plant  variation  within  the  flight  envelope  and  the  eflfects  of  dam¬ 
age  to  aircraft  control  surfaces.  Mathematical  models  of  control  surface  failure  effects  on  aircraft 
dynamics  are  derived  and  used  to  modify  an  existing  small  perturbation  model  of  the  Lambda. 
The  QFT  technique  is  applied  to  design  a  control  system  utilizing  aircraft  pitch  rate,  roll  rate  and 
sideslip  angle  as  feedback  variables.  The  inherent  cross-coupling  rejection  qualities  of  QFT  and  an 
aileron-rudder  interconnect  are  utilized  to  design  a  control  system  which  results  in  a  coordinated 
flight.  An  outer-loop  autopilot  is  then  designed  around  the  QFT  controller  to  further  assist  turn 
coordination.  Sensor  noise  effects  on  aircraft  states  are  also  analyzed. 
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DESIGN  OF  A  FLIGHT  CONTROLLER  FOR 
AN  UNMANNED  RESEARCH  VEHICLE  WITH 
CONTROL  SURFACE  FAILURES  USING 
QUANTITATIVE  FEEDBACK  THEORY 

/.  Introduction 

1.1  Background 

The  Lambda  Unmanred  Research  Vehicle  (UBV)  is  an  asset  of  the  Control  Systems 
Development  Branch  of  Wright  Laboratory  (WL/FIGS),  Wright-Patterson  AFB,  Ohio.  WL/FIGS 
developed  Lambda  to  test  new  flight  control  concepts  and  devices.  Because  tests  on  manned  vehicles 
are  expensive,  WL/FIGS  designed  Lambda  to  offer  a  low  cost,  safe  alternative  to  such  flights. 
When  WL/FIGS  flight  tests  a  new  flight  system,  pilots  located  on  the  ground  send  commands  to 
the  Lambda  via  radio  transmitter.  The  output  of  the  sensors  on-board  the  vehicle  are  transmitted, 
through  telemetry,  to  a  ground  station  in  order  to  be  recorded  for  subsequent  analysis.  New  flight 
control  methods  are  tested  on  Lambda  with  the  hope  that,  if  successful,  these  same  robust  control 
systems  may  be  further  developed  and  implemented  on  operational  Air  Force  aircraft. 

Several  successful  aircraft  flight  compensators  (analog  control  system)  and  flight  controllers 
(digital  control  system),  both  denoted  by  G,  have  been  designed  for  the  Lambda  using  a  feedback 
design  technique  called  Quantitative  Feedback  Theory  (QFT).  This  design  method  guarantees 
desired  plant  response  in  the  presence  of  structured  plsmt  parameter  uncertainties.  QFT  has  been 
used  to  achieve  numerous  robust  aircraft  flight  control  system  designs.  Most  designs  using  QFT 
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have  focused  on  structured  uncertainties  associated  with  varying  flight  conditions.  WL/FIGS  has 
investigated  and  is  investigating  the  feasibility  of  utilizing  a  single  QFT  designed  G  that  maintsdns 
the  aircraft  response  within  specified  limits  with  and  without  control  surface  failures.  Since  this 
thesis  involves  a  digital  control  system,  G  is  henceforth  referred  to  as  the  controller. 

A  flight  controller  that  guarantees  performance  under  degraded  conditions  such  as  control 
surface  failures  is  of  great  benefit.  Such  a  controller  can  greatly  reduce  pilot  effort,  increase  aircraft 
stability  margins,  and  aid  in  successful  aircraft  recovery  if  implemented  on  manned  aircraft.  Imple¬ 
menting  a  controller  with  these  qualities  on  a  test  vehicle  is  the  first  step  in  assessing  its  feasibility 
as  a  viable  controller. 

1.2  Problem 

The  objective  of  this  thesis  is  to  design  a  digital  flight  controller  to  be  implemented  on  the 
Lambda’s  microprocessor-based  digital  flight  control  system.  The  controller  is  to  maintain  Lambda 
performance  in  the  yaw,  pitch  and  roll  channels  with  and  without  control  surface  failures.  Control 
surfaces  considered  are  elevators,  ailerons  and  rudders.  Dr.  Meir  Pachter  has  shown  that  control 
surface  failures  alter  sdrcraft  stability  derivatives  in  the  same  manner  as  varying  flight  conditions  [1]. 
This  relation  is  essential  to  allow  the  application  of  QFT  design  methods  to  control  surface  failures. 
The  design  addresses  failures  at  all  points  of  the  flight  envelope  for  the  Lambda.  Parameters 
defining  the  flight  envelope  consist  of  aircraft  center  of  gravity  location,  airspeed,  altitude,  and 
weight. 

A  reasonable  set  of  failure  modes  are  determined.  WL/FIGS  wants  to  determine  which  control 
surfaces  can  fail  and  still  allow  the  controller  to  meet  performance  specifications.  Control  surface 
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failures  are  modeled  as  a  loss  of  surface  area.  This  type  of  failure  is  more  difficult  to  handle  than 
an  actuator  failure  which  simply  freezes  the  control  surface. 

It  is  desired  to  design  a  controller  which  not  only  adequately  controls  the  aircraft  in  the 
presence  of  control  surface  failures,  but  also  maintains  a  coordinated  ffight.  Coordinated  flight  is 
defined  as  flight  which  minimizes  sideslip  angle.  An  outer-loop  autopilot  is  designed  around  the 
QFT  controller  to  perform  the  coordination  function. 

1.3  Summary  of  Current  Knowledge 

Most  mrcraft  require  a  robust  flight  control  system  for  proper  flight  performance.  Since  the 
acquisition  of  the  Lambda  aircraft  in  1986,  WL/FIGS  has  sponsored  several  theses  which  have 
designed  flight  control  systems  for  the  Lambda.  Capt.  David  Wheaton  wrote  a  Masters  thesis 
which  designs  a  flight  control  system  for  the  Lambda  using  digital  Quantitative  Feedback  Theory 
[2].  This  thesis  shows  that  Quantitative  Feedback  Theory  (QFT)  is  a  viable  design  method  for 
aircraft  flight  control  systems.  The  model  Wheaton  used  in  his  thesis,  however,  was  developed  by 
WL/FIGS  using  only  physical  dimensions  of  the  aircraft.  This  model  does  not  reflect  actual  flight 
data  too  accmately.  Lt.  Gerald  A.  Swift,  in  his  Masters  thesis  work,  derived  a  more  accurate 
Lambda  model  using  actual  flight  data  [3].  Capt.  Donald  J.  Lacey  used  this  improved  model  in  his 
distal  QFT  design  of  a  flight  controller  [4].  This  thesis  expands  the  work  of  Lacey  by  considering 
control  surface  failures  as  an  additional  source  of  uncertainty. 

Quantitative  Feedback  Theory  is  a  robust  multi- variable  feedback  control-system  design  tech¬ 
nique  which  guarantees  system  spedflcations  can  be  met  even  with  structured  plant  parameter  un¬ 
certainties  [5].  Dr.  Isaac  Horowitz,  formerly  of  the  Weitzman  Institute  of  Science,  Rehovet,  Israel, 
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developed  QFT  theory  over  the  past  thirty  years.  He  first  formalized  procedures  in  1979.  During 
this  period,  numerous  methods  of  feedback  design  for  multiple-input  multiple-output  (MIMO)  sys¬ 
tems  had  been  investigated.  None,  however,  addressed  the  topic  of  structured  plant  uncertainty, 
which  is  most  certainly  present  in  most  real-world  control  systems.  While  QFT  has  proven  itself 
as  a  valid  design  technique,  perhaps  its  greatest  asset  is  the  elimination  of  much  of  the  guesswork 
associated  with  other  methods.  QFT,  in  fact,  "consists  of  a  succession  of  direct  (no  iterations  nec¬ 
essary)  single-loop  design  steps”  [5].  These  characteristics  maike  QFT  a  natural  choice  for  mrcraft 
flight  control  system  design. 

Many  QFT  designs  have  been  performed  on  aircraft  other  than  the  Lambda.  This  shows 
QFT  to  be  a  flexible  design  process.  Following  is  a  partial  list  of  AFIT  theses  which  used  QFT  in 
the  design  of  flight  control  systems: 

•  Adams,  J.M.  -  Digital  QFT  Design  for  the  AFTI/F-16^  (1988) 

•  Hamilton,  S.W.  -  Digital  QFT  Design  for  an  Unmanned  Research  Vehicle,  (1987) 

•  Kobylarz,  T.S.  -  Flight  Controller  Design  with  Nonlinear  Aerodynamics,  Large  Parameter 
Uncertainty,  arid  Pilot  Compensation,  (1988) 

•  Miller,  R.B.  -  Multiple-Input  Multiple-Output  Control  System  Design  for  the  YF-16  Using 
Nonlinear  QFT  and  Pilot  Compensation,  (1990) 

•  Neumann,  K.N.  -  A  Digital  QFT  Design  for  the  Control  Reconfigurable  Combat  Aircraft,  (1988) 

•  Ott,  P.T.  -  URV  Reconfiguration  using  Digital  QFT,  (1988) 

•  Rasmussen,  S.J.  -  Application  of  non-linear  QFT  to  Flight  Control  Design  for  High  Angle  of 
Attack  Maneuvers  with  Thrust  Vectoring,  (1991) 
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•  Russdl,  H.H.  -  Analog  QFT  Design  for  the  KC-ISS,  (1984) 


1.4  Assumptions 

It  is  assumed  that  a  linear  time-invariant  (LTI)  minimum-phase  (m.p.)  model  for  each  flight 
condition  accurately  reflects  Lambda  performance.  This  limits  aircraft  maneuvers  to  low  angles-of- 
attack.  The  Lambda  is  designed  to  fly  at  low  angles-of-attack  so  this  assumption  is  valid.  In  using 
an  LTI  model,  the  following  conditions  are  assumed  [6]: 

•  Time  spans  are  short  (30  -  60  seconds)  and  perturbations  small. 

•  The  aircraft  is  assumed  to  be  a  rigid  body  (mass  and  mass  distribution  are  constant). 

•  The  reference  coordinate  system  used  is  a  rotating  system  centered  at  the  aircraft  center  of 
mass. 

•  In  the  horizontal  plane,  the  aircraft  is  symmetrical  about  the  center  of  the  fuselage  (no  cross¬ 
coupling  of  the  longitudinal  and  lateral  directional  equations  of  motion). 

This  design  also  assumes  the  following  properties  of  the  Lambda: 

•  All  inputs  (control  surface  deflections)  are  measurable. 

•  Aircraft  orientation,  rate  and  sideslip  an^e  signals  are  available. 

•  A  digital  sampling  rate  of  50  Hz  is  used  for  the  Lambda  flight  control  system. 

1.5  Scope 

The  objective  of  this  thesis  is  to  develop  controller  difference  equations  which  satisfactorily 
control  Lambda  over  the  entire  flight  envelope,  allowing  for  possible  partial  control  surface  failures. 
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The  design  ends  when  the  computer  simulations  show  Lambda,  with  partial  control  surface  failures, 
responding  within  specifications,  to  control  surface  inputs. 


1.6  Standards 

The  controller  is  designed  to  meet  frequency  domain  tracking  specifications  set  by  WL/FIGS. 
Tracking  specifications  are  given  in  the  frequency  dommn  and  are  shown  in  Tkble  1.1,  where 
is  the  upper  bound  of  the  closed-loop  system’s  transfer  function  and  T/  is  the  lower  bound.  The 
controller  is  required  to  maintain  a  45  degree  phase  margin  angle  for  all  responses. 

It  may  be  difficult  or  impossible  to  meet  the  given  specifications  for  certsun  control  surface 
failures.  In  these  instances,  the  controller  will  guarantee  aircraft  plant  stability  only. 


Table  1.1  Upper  and  Lower  Tracking  Bounds 


Channel 

Tu 

f,  J 

pitch 

106 

^36  II 

»*+St+10O 

«*4-21«>fl44<+330  II 

roll 

II 

yaw 

T 

Ts  n 

,»+7.4M>+l9.St+lS  H 

1.7  Approach/Methodology  I 

Any  control  system  design  ^-^gins  with  the  generation  of  accurate  mathematical  models.  As  I 

part  of  his  thesis  work.  Swift  wrote  several  Matlab  macros  to  generate  plant  models  for  various  I 

flight  conditions.  These  macros  will  be  used  to  generate  plant  data  for  this  design.  Swift’s  macros  y 

were  not  written  to  generate  plants  with  control  surface  failures.  Therefore,  the  macros  are  modified  I 

to  account  for  the  change  in  stability  derivatives  caused  by  the  control  surface  failures.  I 
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The  lon^tudinal  and  lateral-mode  models  of  the  aircraft  are  decoupled  when  no  failures  are 
present,  allowing  the  longitudinal  and  lateral  designs  to  be  performed  separately.  Failure  modeling 
will  attempt  to  maintain  a  decoupled  design.  The  longitudinal  design  reduces  to  a  single-input, 
sin^e-output  (SISO)  design  while  the  lateral  design  becomes  a  two-by-two  multiple-input,  multiple- 
output  (MIMO)  design.  Elevator  and  pitch  rate  are  the  inputs  and  outputs,  respectively,  for  the 
longitudinal  design.  Aileron  and  rudder  are  the  lateral  inputs  and  roll  rate  and  sideslip  angle  are 
the  lateral  outputs.  QFT  design  procedures  vary  slightly  for  these  two  cases. 

The  procedure  for  the  SISO  design  case  is  simpler  than  the  MIMO  case.  The  MIMO  design 
is  actually  a  series  of  SISO  deugns,  one  SISO  design  for  each  input  to  the  system  [7].  The  steps 
for  the  entire  design  process  are  listed  below: 

•  Determine  how  failures  will  affect  stability  derivatives 

•  Alter  Swift  macro  appropriately 

•  Determine  desired  flight  conditions  (including  control  surface  failures) 

•  Generate  aircraft  plants 

•  Form  pseudo-continuous  time  equivalent  transfer  function 

•  Generate  frequency  response  data 

•  Plot  templates 

•  Alter  flight  conditions  as  necessary  to  fill  out  template  (repeat  above  4  steps  until  satisfied) 
At  this  point,  the  SISO  and  MIMO  design  procedures  begin  to  differ.  The  SISO  procedure  is  listed 
first  [8]: 

•  Choose  nominal  plant 

•  Form  stability  and  performance  bounds 
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•  Determine  which  plants  can  meet  performance  bounds 

•  Shape  nominal  loop  (design  compensator) 

•  Synthesize  prefilter 

•  'Dansform  compensator  and  prefilter  to  the  Z-plane 

•  Simulate  using  designed  compensator  and  prefilter 
MIMO  design  procedure  [7]: 


•  Choose  wmghting  matrix  W 

•  Determine  the  plant  transfer  function  matrix,  P(s) 


•  Invert  ?(&),  viz., 


p-'(.)  = 


rtiW  «»(») 


Calculate  Q(s), 


Q(s)  = 


fii(»)  Fu(») 


•  Plot  Templates  for  one  input-output  pair 

■  Choose  a  nominal  plant 

•  Form  stability  and  performance  bounds 

•  Shape  nominal  loop  (design  first  element  of  compensator) 
e  Synthesize  first  element  of  the  prefilter 

•  Shape  the  remaining  loops 

•  Form  the  remaining  compensator  and  prefilter  elements 

•  IVansform  compensator  and  prefilter  to  the  Z-plane 

■  Simulate  using  designed  compensator  and  prefilter 
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1.8  Materials  and  Equipment 


Computer  analysis  for  this  thesis  is  performed  on  the  San  SPARCstation  2  workstations 
available  at  both  the  Air  Force  Institute  of  Techn<dogy  and  Wright  Laboratories.  Software  packages 
used  are  Matlab  with  Simulink  and  the  AFIT  MIMO/QFT  CAD  package. 

1.9  Other  Support 

Dr  Meir  Pachter,  Dr  Constantine  Houpis,  and  Capt  Stuart  Sheldon  have  all  served  on  the 
thens  committe.  Dr  Pachter  and  Dr  Houpis  are  committee  co-chairmen.  Mr  Richard  Sating  assisted 
in  using  the  AFIT  MIMO/QFT  CAD  package. 
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II.  The  Lambda  Aircraft 


t.l  The  Aircraft 

The  Ljtmbda  URV,  shown  in  Fig.  2.1,  is  a  small  remotely-piloted  airplane  with  a  14  foot 
wing  span  and  weight  of  approximately  200  pounds.  It  uses  a  pusher  propeller  and  a  conventional 
horizontal  tail.  The  horizontal  tail  consists  of  a  horizontal  stabilizer  and  a  split  elevator.  The 
vertical  tail  is  located  on  either  end  of  the  horizontal  tail,  and  consists  of  a  vertical  stabilizer  and 
rudder.  The  wings  are  slightly  tapered;  each  has  three  trailing,  movable  control  surfaces.  See 
Table  2.1  for  Lambda’s  characteristics. 

Lambda  has  ten  control  surfaces.  They  are: 

•  Left  Elevator 

•  Right  Elevator 

•  Left  Rudder 

•  Right  Rudder 

•  Left  Aileron 

•  Right  Aileron 

•  Left  Outer  Flap 

•  Right  Outer  Flap 

a  Left  Inner  Flap 

•  Right  Outer  Flap 
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Table  2.1  Lambda’s  Physical  Characteristics 


Dimensioiis  Wing  Span .  14  ft 

Wing  Area .  19  sq  ft 

Length .  9.6  sq  ft 

Height  With  Landing  Gear .  3  ft 

_ Propeller  Diameter . . 2.3  ft _ 

Weights  Maximum  Fuel .  14  lb 

Maximum  Payload .  15  lb 

_ Maximum  Flight  Weight . . 210  lb _ 

Performance  Maximum  Level  Speed  at  Sea  Level.  115  mph  (100  knots) 

Stall  Speed .  63  mph  (55  knots) 

_ Stall  Speed  With  Flaps . . 52  mph  (45  knots) 

Engine  Power .  18  hp 

_ Type .  . 2  cycle,  2  cylinder 

Control  Limits  Elevator  Deflection  Limit .  ±15  deg 

Rudder  Deflection  Limit .  ±25  deg 

Flap  Deflection  Limit .  20  deg  down 

_ Aileron  Deflection  Limit . . ±15  deg _ 


Although  each  control  surface  can  operate  independently,  flight  data  used  to  create  the  mathemat¬ 
ical  model  upon  which  this  design  is  based  are  obtained  with  all  four  of  the  flap  control  surfaces 
operating  together,  as  are  the  rudders.  The  split  elevators  were  deflected  together  while  the  ailerons 
were  operated  diflTerentially.  Flaps  are  not  used  as  a  control  surface.  Thus,  the  region  of  flight  in¬ 
vestigated  in  this  thesis  is  during  cruise  flight  conditions  since  flaps  are  not  involved  as  control 
surfaces  at  this  point  of  flight. 


The  Model 

Lt  Gerald  Swift  developed  a  small  perturbations  model  of  the  Lambda  based  on  data  gathered 
during  a  series  of  flight  tests  [3].  The  model  first  calculates  the  standard  aircraft  stability  derivatives 
for  the  Lambda  based  on  the  desired  flight  condition.  The  derivatives  are  then  placed  into  a  state- 
space  representation.  The  model  is  generated  with  a  Matlab  macro  written  by  Swift.  Capt  Steve 
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Rasmussen  of  Wright  Laboratories  later  modified  some  derivatives  in  the  model  to  better  fit  recent 
flight  tests.  The  model  describes  the  aircraft  dynamics  with  eight  states,  four  states  for  both  the 
longitudinal  and  lateral  channels.  The  ori^nal  model  is  derived  assuming  that  the  longitudinal  and 
lateral  channels  are  decoupled.  The  longitudinal  states  are  u,a,q,  and  9.  Lateral  states  are  P,p,4> 
and  r. 

Swift’s  macro  generates  an  aircraft  model  for  an  undamaged  aircraft.  When  control  surface 
damage  occurs,  the  model  changes  according  to  the  degree  of  damage.  The  mathematical  derivation 
of  these  changes  is  presented  in  Chapter  III.  The  changes  are  incorporated  into  Swift’s  macro.  The 
modified  macro  which  calculates  damaged  aircraft  models  is  presented  in  Appendix  A. 

The  state-space  representation  generated  by  the  macro  is  of  the  following  general  form: 


x(t)  =  A(t)x(t)  4- B(t)u(«) 

(2.1) 

y{t)  =  Cx(t) 

(2.2; 

where  x(t)  is  an  8il  vector  of  aircraft  states,  u(/)  is  a  3x1  vector  of  inputs.  A(t)  and  B(t)  are  8x8 

and  8x3  matrices  respectively.  Because  the  output  of  the  system,  y(t),  is  simply  the  aircraft  states, 

C(t)  is  an  8x8  identity  matrix. 

Assuming  zero  initial  conditions  (xo  =  0)  and  trimmed  flight  (x'o  = 

of  the  plant  is  the  perturbation  model: 

0),  the  Laplace  transform 

sX{s)  =  A(s)X(s) -1- B(s)U(s) 

(2.3) 

Y(s)  =  C(s)X(5) 

(2.4) 
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where: 


u 
a 
9 

e 

/? 

P 

4> 

r 

2.S  Actuators  and  Sensors 

Three  actuator  models  are  used  in  the  mathematical  model  of  the  Lambda;  one  each  for  the 
elevator,  aileron  and  rudder.  Only  one  actuator  model  is  necessary  for  each  pair  of  control  surfaces 
since  each  pair  operates  in  tandem.  Models  were  constructed  by  Wright  Laboratories  based  on 
actual  flight  data.  The  elevator  model  is: 


136.2 

s  +  10.81  ±  j7.74 


(2.8) 


The  aileron  model  is: 


216.84 

s  +  11.07  ±  jTO.09 


(2.9) 
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The  rudder  model  is: 


3.906 
s  +  4.4372 


(2.10) 


Note  that  the  D.C.  gain  ^  1 

Sensors  used  to  measure  aircraft  states  are  not  modeled  explicitly  in  Swift’s  original  macro. 
Therefore,  the  dynamics  of  all  the  sensors  are  included  in  the  model. 
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III.  Control  Surface  Failure  Modeling 


When  damage  to  control  surfaces  occurs,  an  aircraft  loses  the  ability  to  generate  mo¬ 
ments  about  its  axes.  The  aircraft  also  experiences  a  change  in  its  dynamic  response  because 
the  loss  of  an  aerodynamic  surface  changes  the  stability  derivatives  of  the  aircraft.  This  chapter 
derives  the  theory  which  describes  the  changes  in  stabiUty  and  control  derivatives  when  control 
surface  damage  occurs.  The  derived  equations  are  then  implemented  in  the  Matlab  macro  used 
to  generate  aircraft  models.  The  original  macro  was  written  by  Lt  Gerald  Swift  in  his  thesis. 
The  original  macro  calculates  all  aircraft  dimensional  and  non-dimensional  stability  derivatives  for 
an  undamaged  aircraft  at  jpven  flight  conditions.  The  ‘Tailed”  stability  derivatives  (denoted  with 
subscript  F)  are  expressed  in  terms  of  the  “unfailed”  stability  derivatives  (denoted  with  subscript 
UF). 


3.1  Longitudinal  Derivatives 

The  Lambda  has  two  sets  of  lon^tudinal  control  surfaces;  the  elevators  and  flaps.  Flaps  are 
used  primarily  for  take-off  and  landing  maneuvers.  Because  this  thesis  investigates  flight  control 
at  cruise  flight  conditions,  only  the  elevator  is  considered. 

The  variable  is  defined  as  the  damage  level  to  the  elevators.  Damage  to  be  considered  is 
in  the  form  of  lost  surface  area.  Such  damage  on  an  aircraft  might  occur  as  the  result  of  enemy 
fire.  Ce  is  the  percentage  of  elevator  surface  area  which  remains  after  being  damaged.  When  Ce  =  1, 
no  damage  has  occurred.  When  Ce  =  0  both  elevators  are  completely  blown  off.  If  elevators  are 
damaged  equadly  no  coupling  to  the  lateral  channel  exists.  Unequal  damage,  though,  creates  a 
lateral  moment.  For  modeling  purposes,  this  analysis  cissumes  only  one  elevator  is  damaged  at  any 
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one  time.  Therefore,  the  ramge  0.5  <  C«  ^  thesis.  Ensuing  chamnel  coupling  is 

discussed  in  more  detail  in  Section  3.3.1. 

The  effects  of  failure  can  be  determined  by  analyzing  the  three  non-dimensionalized  longitu¬ 
dinal  equations  of  motion  as  given  in  Blakelock  [Eq  (1-59)]. 

<  W*  ■  (3‘) 

-(C,.u)  +  [(^  -  ^c,ja]  -  C,.a  +  ((-^  -  ^C,,)*  -  C.(Mne)«l  =  (3.2) 

-(C„.ti)  +  (-^C„,a  -  C„.o)  +  (i*  -  ^C„J)  =  (3,3) 

The  Ci-  terms  are  the  non-dimensional  stability  derivatives.  These  derivatives  describe  the  longi¬ 
tudinal  dynamics  of  the  aircraft.  They  will  change  when  elevator  damage  occurs. 

Some  simplifying  assumptions  are  made.  Fuselage  effects  are  assumed  negligible.  Also,  effects 
on  the  derivatives  in  Eq.  (3.1)  can  be  assumed  negli(pble  as  changes  in  elevator  surface  area  create 
negligible  forces  in  the  x  direction.  This  leaves  only  derivatives  in  the  Z-equation  [Eq.  (3.2)]  and 
the  M-equation  [Eq.  (3.3)]  to  be  analyzed.  The  effect  on  C,„  can  also  be  neglected  because  this 
term  is  related  to  gravity  only. 

3.1.1  M  Equation  Analysis.  Failure  effects  on  several  derivatives  in  Eq.  (3.3)  are  assumed 
negligible:  Cm,  is  due  mostly  to  lift  on  the  wing.  Elevator  dam^e  has  negligible  effect  on  this 
term.  Because  the  final  control  system  will  still  have  a  pilot  in  the  loop,  only  short  term  dynamics 
are  a  concern.  Partials  of  u  are  only  present  in  phugoid  dynamics.  For  this  reason  also,  failure 
effects  on  Cm,  can  be  neglected.  This  assumption  is  based  on  the  fact  that  no  spillover  of  phugoid 
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poles  occurs  when  loops  are  closed.  Final  simulations  will  verify  that  no  spillover  of  short-period 
and  phugoid  dynamics  exists. 

accounts  for  downwash  effects.  It  also  changes  negligibly  when  elevator  damage  occurs. 
Unsteady  dynamics  associated  with  d  are  neglected. 

The  above  discussion  leaves  three  short-period  derivatives  which  must  be  modeled  for  failures, 
C'm,  and  Cm,.-  The  change  in  is  derived  first.  Two  forces  cause  moments  about  the 
center  of  gravity.  These  moments  are  caused  by  the  lift  on  the  wing  and  the  lift  on  the  tail.  Each 
moment  can  be  expressed  as: 


Mwin§  •“  r  S^d^^XcgOt 

Ig 

(3.4) 

Mtaii  =  -  fc)a 

*9 

(3.5) 

Variables  and  AXe* 

have  opposite  signs  for  the  Lambda  because  they 

are  measured  aft  and 

forward  of  the  center  of  gravity,  respectively.  The  total  moment  about  the  center  of  gravity  is 

therefore: 

—  Mwing  "i"  ^tail 

Mtotai  =  y-[5,„au,AA’c,  +  vStaM^  -  k)]a 

(3.6) 

Define  Maup  as: 

=  Y[Sy,a^^Xcg  +  T?5,a,L(l  -  k)] 

(3.7) 

where  is  a  dimensional  stability  derivative  for  the  undamaged  aircraft. 


If  damage  occurs  to  the  elevator,  less  moment  is  produced  by  the  tail/elevator  combination. 
The  change  in  moment  is  proportional  to  the  effective  area  of  the  horizontal  tail.  Therefore  £q.  (3.5) 


becomes: 

*9 


(3.8) 


where 


.  (C«  -  1)5.  +  St 

St 


(3.9) 


C«  is  the  ratio  of  available  tail  surface  area  to  undamaged  tail  surface  area.  When  Ce  =  1  i^his  ratio 
is  1,  thus  restdting  in  a  value  for  Mtan  being  the  same  as  for  an  undamaged  aircra. 

With  a  damaged  elevator,  the  dimensional  stability  derivative  {Ma^)  becomes 


+  (^ftfSttiilt{l  —  Af)]  (3.10) 

Because  C«  is  less  than  1  when  damage  has  occurred,  it  is  seen  that  the  magnitude  of  the  moment 
about  the  center  of  gravity  is  decreased  when  an  elevator  loses  surface  area,  as  expected. 

The  effect  on  can  be  expressed  in  terms  of  the  change  in  the  dimensional  stability 
derivative  AAf^  where 

AMo  =  —  Maf  (3.11) 

Substituting  Eq.  (3.7)  and  (3.10)  mi-o  (3.11)  yields 


Ai»/„  =  f[(l-C>W.(l-*)] 

*9 


(3.12) 
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Td  non-dimensionalize  this  chaoge,  divide  Eq.  (3.12)  by  £q.  (3.10)  to  yield 


AM,  ^  (1-C:)i?5,a,/«(1-A:) 

Af„  +  CnS,aMl  -  k)  ^ 

This  ratio  represents  the  fraction  of  which  has  been  lost  due  to  damage.  Therefore,  the  first 

damaged  dimensionless  stability  derivative  is 


(3.14) 


The  terms  are  subtracted  because  elevator  damage  causes  a  larger  moment. 

Figure  3.1  shows  how  varies  with  damage  fraction  Ce  for  a  representative  flight  condition. 

Figure  3.1  Mo  vs.  C, 


The  two  remaining  stability  derivatives  of  Eq.  (3.3)  are  relatively  simple  to  derive.  The  first, 
is  related  to  the  tail  surface  area.  depicts  damping  in  pitch  motion.  Because  pitch 
damping  is  performed  by  the  tail,  the  derivative  is  directly  related  to  efiective  tail  surface  area 
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therafdra  a  direct  proportional  function  of  . 

=  CCn^,  (3.15) 

has  a  similar  simple  relation.  In  this  case  though,  the  change  in  the  derivative  is 
directly  related  to  the  portion  of  total  elevator  area  which  has  been  lost,  not  to  the  effective  area 
of  the  tail. 

(3.16) 

S.l.t  Z  Equation  Anaiyoia.  Forces  in  the  z-axis  direction  are  described  in  Ek}.  (3.2)  which 
is  repeated  here: 

-  (C,.«)  +  ((^  -  5J7C,»  -  C,.a  +  ((-^  -  ^C,.)«  -  C,(«ne)»l  =  C,J.  (3,17) 

For  the  same  reasons  given  for  the  M  equation,  failure  effects  on  the  partials  of  u  and  d  are 
neglected.  C.  is  again  neglifpble  as  it  is  related  to  gravity  only. 

The  Cj.  term  is  similar  to  the  Cm.  term.  Force  in  the  z-axis  is  generated  by  lift  on  the  wing 
and  the  tail.  Both  forces  are  negative  since  they  act  in  the  negative  z  direction. 

^win$  ~  (3.18) 

Zt,i,  =  --^175,0,(1-*)  (3.19) 

m 

The  sum  of  these  two  forces  yields  Zavp 

Zav,  =  -  *)]  (3-20) 
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Elevator  damage  affects  the  force  caused  by  the  tail  only.  The  force  is  proportional  to  the 


^ective  area  of  the  horizontal  tail.  Therefore  is  shown  as: 


Zap  =  -  — (S*0«  +  CeVStOtil  -  *)] 
m 


(3.21) 


The  difference  in  Za  for  the  failed  and  unfailed  cases  is  AZa  —  Zaup  Zap- 


AZa  =  -Mi~OvSMi~k) 

m 


(3.22) 


The  above  relation  is  divided  by  Eq.  (3.21)  to  yield: 


AZg  (l-QqWl-fc) 
Za  Sp,ap,  +  C«»?5,o,(l  -  k) 


(3.23) 


This  ratio  is  the  fraction  of  C,.  which  has  been  lost  to  damage. 


C...  =  c..„  + 


(3.24) 


When  the  aircraft  experiences  a  pitching  motion,  changes  in  lift  are  generated  primarily  by 
the  tail.  is  therefore  proportional  to  the  effective  horizontal  tail  area. 


(3.25) 


The  tail  effectiveness  is  directly  related  to  available  elevator  surface  area. 


(3.26) 


S.2  Lateral  Derivatives 


Similar  to  the  longitudinal  channel,  the  Lambda’s  lateral  channel  has  two  sets  of  control 
surfaces.  Lateral  control  surfaces  are  the  rudders  and  ailerons.  Both  sets  are  used  during  routine 
flight  to  perform  aircraft  maneuvers.  Both  sets  are  considered  in  damage  analysis. 

Two  damage  variables  are  defined  for  the  lateral  channel,  one  for  each  set  of  control  surfaces. 
Damage  entails  a  loss  of  surface  area  and  is  proportional  to  the  latter.  Damage  to  the  ailerons  is 
denoted  by  Rudder  damage  is  similarly  Cr*  Lateral  control  surface  damage  variables  account  for 
damage  to  both  surfaces  in  each  pair.  For  example,  if  both  rudders  have  lost  25%  of  thar  surface 
area,  the  total  lost  is  50%.  This  means  (r  is  0.5.  Lateral  damage  can  be  represented  as  such  because 
damage  to  a  right  or  left  lateral  surface  has  the  same  net  effect  on  aircraft  dynamics.  No  coupling 
of  lateral  damage  to  the  longitudinal  channel  exists  (see  Section  3.3.1). 

Failure  analysis  begins  with  the  three  lateral  equations  of  motion.  Written  in  Laplace  trans¬ 
form  form,  from  Blakelock  [Eq.  (3-17)]: 

-  C,.#+  (^  -  -  C,.0+  -  C,.0  =  C,..S.  +  C,.A  (3.27) 

-  w*  - 

”  ~  ~  ~  (3.29) 

It  is  seen  that  stability  derivatives  exist  for  the  V*,  P,  r  and  /d  states.  These  derivatives  describe 
the  lateral  dynamics  of  the  aircraft.  These  will  change  when  aileron  and/or  rudder  damage  occurs. 

Like  the  longitudinal  channel,  fuselage  effects  on  the  control  derivatives  are  assumed  negligible. 
C'/.n*  and  Ci,„^  denote  the  change  in  the  respective  moments  due  to  a  change  in  aircraft  orientation. 
These  four  derivatives  are  always  zero.  They  are  not  considered  in  this  analysis. 
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S.t.l  Y  Equation  Analysis.  The  Y  equation  is  repeated  here: 


-  -  C..*  +  (3.30) 

The  failed  derivatives  for  p,  r  and  /3  are  directly  proportional  to  the  unfailed  derivatives. 
These  derivatives  measure  the  force  in  the  y  direction  created  by  the  respective  state.  Only  the 
^ective  area  of  the  vertical  fins  creates  a  force  in  the  y  direction.  Therefore: 


~  Cr^twVP 
^WtP  ~  Cr^SrVP 

^t$p  ~  Cr^tK/P 


(3.31) 

(3.32) 

(3.33) 


where 

_  {(r  -  1)5,  +  5; 

is  the  effective  area  of  the  vertical  fin  with  rudder  damage. 

The  two  derivatives  for  the  control  surfaces  themselves  are  proportional  to  the  area  lost  on 
each  control  surface.  As  each  of  these  derivatives  is  a  measure  of  control  surface  effectiveness,  the 
effectiveness  is  reduced  proportionally  to  the  available  control  surface  area. 


(3.35) 

(3.36) 

For  the  aircraft  model  being  used,  =  0  for  all  cases. 


^ytrP  ~  CrCytrVP 
^ytmP  ~  Ca^t§mUP 
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S.i.S  L  Equation  Analysis. 


The  L  moment  equation  [£q.  (3.28)}  is  repeated  here: 


J2_ 

Sqb 


6  -  —Cl  6  -  —ib  - 
^  2U  Sqb^ 


b 

2V 


+  Ci^^Sr 


(3.37) 


The  C{,  derivative  is  proportional  to  lateral  control  surface  failure.  This  derivative  is  unique 
in  that  it  can  be  altered  by  failures  in  either  the  aileron  or  rudder.  This  results  from  the  fact  that 
both  the  wings  and  vertical  fins  damp  rolling  motion.  The  proportion  factor  Car  is  the  effective 
area  of  the  wings  and  vertical  fins. 


(Ca-l)5»  +  (Cr-l)5,  +  5^  +  5y 
S*  +  Sf 


(3.38) 


The  failed  value  of  C|,  is  therefore 


~  ^ar^lfvr 


(3.39) 


The  other  two  rate  derivatives  are  both  proportional  to  the  effective  tail  area.  Only  the 
vertical  tails  provide  damping  due  to  r  or  movement. 


a,  =  (3.40) 

Clap  =  CCia^p  (3.41) 


The  control  surface  derivatives  are  directly  related  to  their  respective  control  surface  areas. 


Clt»p  — 

(3.42) 

CUtF  ~  ^aCla^fjp 

(3.43) 
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S.2,S  N  Equation  Analysis.  The  N  moment  equation  [Eq.  (3.29)]  is  repeated  here: 


(3.44) 


All  of  the  N  equation  rate  derivatives  are  proportional  to  the  effective  vertical  fin  area  Cr- 
The  vertical  fin  is  the  only  surface  which  has  an  effect  on  moments  about  the  2:-axis. 


=  Cr^nrUF 

^n,p  — 

^n0p  —  Cr^npuF 


(3.45) 

(3.46) 

(3.47) 


The  control  surface  derivatives  are  directly  related  to  their  respective  control  surface  areas. 


^nf»r  —  Cr^i 


nttvr 


^nirP  ~  Ca^'i 


ntrUF 


(3.48) 

(3.49) 


3-11 


S.S  Disturbance  Modeling 


The  lateral  and  lonfptudinal  models  both  have  an  external  disturbance  input.  The  origins  of 
these  disturbances  are  described. 

As  shown  by  Eq.  (2.3),  the  equations  of  motion  for  both  the  longitudinal  and  lateral  channels 
are  written  in  state-space  form  as 


5X(s)  =  A(s)X(s)  -I-  B(s)U(s)  (3.50) 

where  X  is  the  state  vector  for  the  channel  in  question.  Because  the  model  is  time  invariant,  the 
A  and  B  matrices  are  not  functions  of  s  and  they  are  hence  so  defined. 

Eq.  (3.50)  is  augmented  to  include  disturbances  as 

sX(s)  =  AX(s)  -I-  BU(s)  -I-  Td{s)  (3.51) 

where  F  is  a  disturbance  input  matrix  and  d  is  the  disturbance  input(8).  In  the  following  sec¬ 
tions,  the  disturbance  input  matrices  F/at  and  F(„,  are  determined  for  the  lateral  and  longitudinal 
channels,  respectively. 

First  consider  damage  to  lateral  control  surfaces.  If  the  rudder  loses  surface  area,  the  rudder 
loses  ability  to  create  a  yawing  moment.  The  effect  on  the  lateral  channel  is  described  in  Sec.  3.2. 
The  configuration  of  the  two  rudders  does  not  allow  it  to  create  a  pitching  moment,  even  when 
one  rudder  is  damaged  more  than  the  other.  Thus  a  rudder  failure  does  not  cause  a  coupling 
disturbance  in  either  the  lateral  or  longitudinal  channels. 
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The  location  of  the  ailerons,  parallel  to  the  y-axis,  allows  the  ailerons  to  create  a  pitching 
moment  if  both  ailerons  do  not  create  the  same  moment  about  the  y-axis.  If  one  aileron  has  less 
surface  area  than  the  other,  unequal  pitching  moments  occur.  However,  aileron  failures  are  not 
modeled  as  coupling  into  the  longitudinal  channel  because  the  moment  caused  by  the  ailerons  is  a 
small  fraction  of  the  moment  caused  by  the  wing.  Tail  moments  already  counter  the  wing  moments. 

Therefore,  only  elevator  failures  are  modeled  as  disturbances. 

S.S.l  Lateral  Disturbance.  Swift’s  models  for  an  undamaged  aircraft  assume  that  the 
longitudinal  and  lateral  modes  are  uncoupled.  This  is  a  valid  approximation  for  the  undamaged 
aircraft.  However,  the  possibility  of  coupling  must  be  considered  for  a  damaged  aircraft.  Loss  of 
surface  area  may  cause  paired  control  surfaces  to  exert  unequal  moments  about  the  other  channel’s 
axes. 

The  effect  of  damage  to  the  elevator  is  considered.  The  left  and  right  elevator  of  the  Lambda 
are  located  on  the  horizontal  tail  trailing  edge.  The  elevators  share  a  common  inside  edge.  The 
elevators  move  together  when  deflected  to  produce  a  moment.  Because  the  elevators  are  symmetric 
about  the  x-axis  of  the  aircraft,  no  lateral  moments  are  produced.  If  one  elevator  operates  with 
less  surface  area,  it  produces  less  moment  about  the  x-axis  than  the  undamaged  elevator.  The  net 
result  is  a  coupling  into  the  lateral  channel  in  the  form  of  a  moment  about  the  x-axis.  This  moment 
is  considered  to  be  significant  and  therefore  is  included  in  the  aircraft  model. 

The  coupling  of  an  elevator  failure  to  the  lateral  mode  is  modeled  as  a  disturbance  to  the 
lateral  states.  By  modeling  as  a  disturbance,  the  designed  QFT  controller  attenuates  the  effect  of 
this  moment  on  the  output,  mmntaining  an  essentially  uncoupled  system. 
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An  elevator  deflection  which  contributes  an  z-axis  moment  does  so  in  much  the  same  manner 
as  an  suleron  deflection.  The  elevator  ” effectiveness”  can  be  introduced  into  Eq.  (3.28).  This 
additional  term  is  the  moment  about  the  z-axis  added  by  an  elevator  deflection.  In  this  case,  6^  is 
the  elevator  deflection  from  horizontal,  not  from  trim.  Note  that  =  ^=:f,  ^  =  p  and  ij;  =  r. 


(3.52) 


The  above  equation  is  placed  into  state  space  form.  This  requires  solving  Eq.  (3.52)  for  p  as 
follows: 


5  +  C,,J,  +  (3.53) 


Equation  (3.54)  is  now  written  in  terms  of  its  dimensional  derivatives. 


Jtz 

P  ~  ^pP  +  ~r~^  +  ^fa^o  +  Ltr^r  +  LfeSg 


(3.55) 


The  L  and  N  equations  are  coupled  equations.  Therefore,  to  place  in  state  space  form,  f  must 
be  expressed  in  terms  of  the  other  states  and  p.  This  equation  is  then  substituted  into  Eq.  (3.55). 
Notice  no  ^nt,  term  is  included.  Elevator  deflection  is  assumed  to  create  a  negligible  yaw  moment. 
Equation  (3.29)  is  now  solved  for  f. 


Sqb^  Sqb^  2t^^"'’’ 

.  ,  Sqb^  SqP  .  A  A 


(3.56) 

(3.57) 

(3.58) 
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This  equation  is  now  written  in  terms  of  dimensional  derivatives: 


r=^P+N,p+  Nrt  +  Np0  +  Ntah  +  Ntrir 
Substituting  £q.  (3.59)  into  Eq.  (3.55)  and  solving  for  p  yields: 


(3.59) 


AAA  A  A  A 

(3.60) 


•3 

where  kj  =  J„//,  and  A  =  1  -  j^. 


Coupling  exists  between  the  L  and  N  moment  equations  when  placed  in  state-space  form. 
This  is  evidenced  by  the  presence  of  both  f  and  p  terms  in  Els.  (3.28)  and  (3.29).  A  derivation 
similar  to  the  one  just  performed  but  beginning  with  the  N  equation  yields  a  disturbance  which 
affects  the  r  state.  Substituting  Eq.  (3.55)  into  Eq.  (3.59)  yields: 


[1  — =  -^[Lpp-\- LfT LpP Lta6a-\- Lfr6r-\' Npp-^NfT -V N N N tr^r  (3.61) 

r  =  ( - - )p-|-( - - )r-f-( - - )^-K - - )da-l-( - - )gr-l-(  ^  )fl» 

(3.62) 


where  Kn  =  Jgz/Iz  and  A  = 


A  valid  approximation  for  £<*  is  now  derived.  The  elevator  and  aileron  have  approximately 
the  same  area.  They  can  therefore  be  assumed  to  create  the  same  force  when  deflected  by  equal 
amounts.  The  moment  created  about  the  x-axis,  however,  is  different  due  to  different  moment 
arms.  The  elevators  are  approximately  |  the  distance  from  the  x-axis  as  the  ailerons.  Returning 
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to  the  non-  'limensional  unfailed  derivatives  allows  the  following  relationship  to  be  derived: 


^U,uF  — 


(3.63) 


and  similarly: 


LieUF  —  ^LiaUF 


(3.64) 


Because  Lu  is  a  measure  of  elevator  "effectiveness'’  in  the  lateral  channel,  Luf  is  a  proportion 
of  LfeUFt  is: 


LseF  =  (1  ~  Ce)I'ftUF 


(3.65) 


which  can  be  written: 


T  _(1-Ce)r 
^'ieF  —  - g - I'tttUF 


(3.66) 


Examination  of  Eq.  (3.60)  and  Eq.  (3.62)  shows  the  contribution  of  an  elevator  deflection  to 
the  lateral  dynamics  as: 


/?■ 

?  =A 


^  I  + 

4> 


(3.67) 


Substituting  Eq.  (3.66)  into  Eq.  (3.67)  yields: 


0a]  \  0]  \  0 

w  P  . 

yi  j. 

»  J  L  J  L  6A 


The  2x1  Fjat  matrix  is  now  written  from  Eq.  (3.68)  as: 


(3.68) 
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r,.i 


0 

<A 

0 


(3.69) 


L  6A  J 

For  disturbance  modeling  purposes,  St  is  set  to  the  maximum  value  of  elevator  deflection. 


S.S.2  Longitudinal  Disturbance.  Elevator  damage  is  modeled  as  a  disturbance  in  the 
longitudinal  channel  also,  dor  to  damage,  the  elevator  is  set  at  a  particular  trim  condition.  When 
damage  occurs,  forces  and  moments  generated  by  the  elevator  trim  setting  are  lost.  Thus,  the 
elevator  moves  to  a  position  different  from  the  trim  position  for  the  linearized  model.  The  change 
in  forces  and  moments  are  therefore  modeled  as  an  elevator  deflection. 

The  moment  generated  by  an  undamaged  elevator  is  given  as: 


Mvf  = 


(3.70) 


where  is  used  to  avoid  confusion  with  the  devator  deflections  of  Ekjs.  (3.1),  (3.2)  and  (3.3). 
When  the  elevator  is  failed,  the  moment  is  a  fraction  of  available  elevator  surface  area: 


Mp  =  CtqStOtSt 


(3.71) 


The  moment  generated  by  an  undamaged  trim  setting  is: 


MuF  =  qStOtStrim 


(3.72) 
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where  Strim  is  referenced  to  the  horizontal  tail  surface.  The  trim  moment  lost  when  elevator  damage 
occurs  is  a  function  of  the  undamaged  elevator  trim  moment: 

Aff  =  (l-C«)gWi«m  (3.73) 

After  damage  occurs,  the  elevator  trim  adjusts  to  make  up  for  the  lost  moments.  Equation  (3.71) 
and  Eq.  (3.73)  are  therefore  set  equal  to  determine  the  required  elevator  deflection. 

=  (1  ~  OgStOtStrim  (3-74) 

K  -  (3.75) 

of  Eq.  (3.75)  is  the  elevator  deflection  required  to  offset  the  lost  trim  moments  caused  by  elevator 
damage.  This  deflection  appears  in  the  lon^tudinal  equations  of  motion  in  addition  to  the  standard 
control  deflection  Sg.  The  effectiveness  of  is  the  same  as  ig.  The  z  and  m— equations  are  rewritten 
with  the  additional  moment  and  force  terms.  The  z— equation  is  neglected  because  =  0. 

=  (3.76) 

-  (C~.“)  +  -  C™.o)  +  +  C^X  (3.77) 

The  additional  moments  and  forces  are  added  because  the  additional  trim  deflection  creates  mo¬ 
ments  and  forces  in  the  same  directions  as  a  normal  dynamic  control  deflection. 
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Similai  to  the  lateral  channel  disturbances,  the  above  equation  needs  to  be  placed  in  state- 


space  form. 

(3.78) 

(3.79) 

Setting  in  terms  of  dimensional  derivatives: 


9  =  +  ( 


=  %=•«  +  + c-o)  + 

59c*  ,  .  Sqc^  .X  .  Sqc^  ^  .x  ,  Sqc ^  .  Sqc^  . 

+  J  Cni«Qi)+  2/  /  ^”***^*  "^  J 


9  =  Af„tt  Afad  +  AfaCt  +  M^q  +  Mt^Sf  +  Af*,Aj 


(3.80) 


Now  E^.  (3.2)  is  used  to  solve  for  d  and  substitute  that  result  into  3.80.  Adding  the  additional 
trim  forces  to  Eq.  (3.2)  yields: 


^(O  -  =  ^«  +  C,.o  +  (^(P  +  ^c.,)g  +  C.(«ne)»  +  C.,.6,  +  C.,j:  (3.82) 


m 


m 


m 


m 


(3.83) 


Now  define  in  non-dimensional  terms: 


[U  —  2ra]d  =  Z^u  -|-  ZaQ  -|-  [17  -j-  Zg]q  —  g{sinQ)6  -I-  Z(^6e  -f-  Zs,^S^  (3.84) 

^  _  _^L_«  4.  n  4.  +  ^<]  _  g(g»'ne)g  f  .  e-x 

“  “  [f7  -  Za]“  [t^  -  Z^]"  +  [tr  -  Za]’  [U  -  Za]  [t^  -  Za]^*  [f7  -  Za]^*  ^  ^  ^ 
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Sabstitnte  Eq.  (3.85)  into  £q.  (3.80)  to  obtain  the  q  state-space  equation: 


(3.86) 

By  extracting  the  terms  in  Eqs.  (3.85)  and  (3.86)  related  to  the  contributions  of  elevator 
damage  to  a  longitudinal  disturbance  can  be  written  in  matrix  form. 


’  " 

'  ttd  ' 

ad 

id 

=  Ap 

Qd 

9d 

+ 

.  »d  . 

.  . 

0 

V-zi 

0 


Substituting  Eq.  (3.75)  into  the  above  yields  ri„«: 


(3.87) 


r  ♦  *1 

’  ' 

Od 

=  Ap 

ad 

id 

9d 

c. 

.  9d  . 

.  . 

0 


Mi, 


'•  “  U-Zi 


_ a. 

Zi. 


^trim 


where 


(3.88) 


r/n,  = 


(1-Ce) 

C. 


0 

JiM. 


V-Zi 

"<.+^ 

0 


(3.89) 


For  simulation  purposes,  the  maximum  expected  value  of  elevator  trim  is  used  as  Strim-  This 
derivation  reveals  that  the  effect  of  lost  trim  moments  can  be  modeled  as  a  disturbance  in  the 
longitudinal  channel.  The  disturbance  is  proportional  to  the  unfailed  trim  setting  and  is  a  function 

of  C*. 
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Figure  3.2  MIMO  QFT  CAD  Package  Block  Diagram 

S.S.S  Plant  and  Disturbance  Matrices.  Figure  (3.2)  shows  the  general  block  diagram  of 
the  model  used  by  the  MIMO  QFT  CAD  Package.  The  derived  failure  models  need  to  be  described 
in  terms  of  the  plant  matrix  P(s)  and  the  disturbance  plant  matrix  Pj(s). 

The  failure  models  have  been  derived  in  the  form: 


X(t)  =  A(t)x(t)  +  B(t)n(t) 

(3.90) 

II 

n 

(3.91) 

Placing  these  equations  into  Laplace  domain  yields: 


sx(s)  =  A(s)x(s)  -(-  B(s)u(5) 

(3.92) 

y(s)  =  Cx(s) 

(3.93) 
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The  statM  ajid  outputs  are  therefore  solved  as: 


x(3)  =  [sl  -  A]-*Bu  +  (si  -  A]-‘rd(s) 


(3.94) 


and 


y(s)  =  C(sl- Aj-^Bu  +  CfsI- A]-‘rd(s) 

=  P(s)u(s)  +  P4(s)d(s)  (3.95) 


where 


P(s)  =  C(sI-A]-*B  (3.96) 

P^(s)  =  C[sI-A]-‘r  (3.97) 

Equation  (3.95)  is  in  the  form  represented  by  the  system  block  diagram  in  Fig.  (3.2).  Thus 
the  lateral  system  matrices  are: 


P(s)  =  C[sI-Af]-^Bf  (3.98) 

Prf(«)  =  C[sl  —  A/-]"'rioi  (3.99) 

where 

dut{s)  =  6,{s)  (3.100) 

and 
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r,.,= 


The  longitudinal  system  matrices  are: 


6A 

0 

«A 


(3.101) 


where 


P(s)  =  C[sI-Af]-‘Bf 


(3.102) 


P,(s)  =  c[si-AH-‘r,„, 


(3.103) 


(3.104) 


r...  =  ^ 

0 


(3.105) 
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IV.  Longitudinal  MISO  Design 


This  chapter  discusses  the  design  of  the  lon^tudinal  controller.  The  design  steps  are  outlined 
in  Chapter  1. 

4.1  Plant  Generation 

4.1.1  Flight  Conditions.  Plants  are  selected  to  adequately  represent/cover  ail  regions  of 
the  Lambda  flight  envelope  and  all  possible  flight  conditions.  It  is  desired  to  include  as  large  a 
variation  in  the  flight  parameters  as  possible.  A  la^e  variation  better  ensures  the  QFT  controller 
will  adequately  control  the  aircraft  at  any  flight  condition. 

Uncertain  flight  parameters  are  chosen  with  consideration  to  the  actual  flight  envelope  of 
the  Lambda  and  to  the  expected  conditions  during  flight.  The  flight  parameters  consist  of  four 
variables: 

•  Center  of  Gravity 

•  Airspeed 

•  Altitude 

•  Weight 

For  the  investigation  of  control  surface  failure  effects  on  the  longitudinal  channel,  elevator  damage 
is  included  as  a  fifth  variable.  It  is  desired  to  reduce  the  range  on  the  variables  when  possible 
without  eliminating  a  flight  condition  which  the  aircraft  might  encounter.  This  reduces  excess 
plant  variation  from  the  QFT  design  which  facilitates  design  of  the  controller. 
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At  air  speeds  less  than  70  knots  the  Lambda  has  unstable  phugoid  poles.  The  time  constant 
of  the  unstable  phugoid  poles  is  very  slow.  These  plants  are  eliminated  from  consideration  as  the 
Lambda  does  not  routinely  operate  at  speeds  less  than  70  knots.  Speeds  less  than  70  knots  require 
use  of  the  flaperons  for  control.  Flaperons  are  not  considered  in  the  design.  The  maximum  speed 
of  the  Lambda  is  100  knots.  To  adequately  represent  upper  airspeed  limits,  a  maximum  speed  of 
110  knots  is  used.  Therefore,  airspeed  is  varied  from  70  knots  to  110  knots. 

Lacey  uses  a  weight  range  of  181  to  215  lbs  for  his  design.  This  range  is  somewhat  narrowed 
for  this  design.  The  weight  of  the  Lambda  without  fuel  is  194  lbs.  With  a  full  tank  at  take-off,  the 
weight  is  208  lbs.  To  encompass  this  range  and  to  allow  additional  components  to  be  loaded  on 
the  Lambda,  a  wdght  range  of  190  to  215  lbs  is  used. 

The  high  altitude  value  of  the  Lambda  flight  envelope  is  10000  ft.  Because  altitude  has 
large  effects  on  the  plant  dynamics,  a  range  of  1000  to  10000  ft  is  considered.  It  is  assumed  that 
high-dynamic  maneuvers  are  not  performed  below  1000  ft. 

Finally,  the  same  center  of  gravity  range  as  used  by  Lacey  is  chosen  for  the  design.  This 
allows  for  c.g.  variation  due  to  fuel  consumption  during  flight.  The  chosen  range  is  45.83”  to 
47.74”  measured  from  the  aircraft  nose. 

Because  the  extent  of  the  effect  of  elevator  failure  on  aircraft  dynamics  is  not  known,  the 
minimum  Ce  Is  chosen  to  be  0.75.  This  corresponds  to  half  of  one  entire  elevator  being  lost.  This 
is  a  substantial  failure.  If  design  efforts  show  that  a  QFT  design  cannot  be  accomplished  for  this 
level  of  a  failure,  Ce  wiU  be  altered  to  reflect  less  damage. 

Each  flight  parameter  is  taken  at  its  extreme  maximum  and  minimum  values.  With  five 
conditions  being  varied,  32  flight  conditions  are  analyzed.  The  number  of  plants  actually  used  in 
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the  CAD  package  is  reduced,  however,  for  computational  ease.  This  is  discussed  in  further  detail 
in  Section  4.4.  The  initial  32  flight  conditions  are  presented  in  Appendix  B. 

4.1.2  Short-Period  Approximation.  To  reduce  computations  and  increase  their  numerical 
accuracy,  the  short-period  approximation  is  considered.  Slow  aircraft  dynamics  are  considered 
negligible.  With  a  pilot  in  the  loop,  pilot  corrections  are  anticipated  that  will  negate  any  slow 
(phugoid)  d3mamics.  Therefore,  if  the  short-period  approximation  varies  from  the  full  longitudinal 
dynamics  only  at  suflSciently  slow  frequencies,  the  short-period  approximation  can  be  considered 
a  valid  approximation.  Figure  4.1  shows  that  the  short-period  approximation  duplicates  the  full 
equations  above  «>  =  1^-  Frequencies  below  are  considered  slow  dynamics.  A  numerical  rule 

Figure  4.1  Frequency  Response  of  Complete  Longitudinal  Eqs.  and  Short-Period  Approximations 


Compfete  Longiiucfnal  Short-Period  Approximation 


of  thumb  for  use  of  the  short-period  approximation  is  that  This  relation  holds  true 

for  all  32  flight  conditions.  Hence,  use  of  the  short- period  approximation  is  considered  valiu  from 
a  theoretical  point  of  view. 

Final  verification  of  the  validity  of  the  short-period  assumption  is  obtained  when  simulations 
are  performed.  Simulations  are  performed  with  the  full  model  of  aircraft  dynamics.  Any  spillover 
of  the  phugoid  into  the  short-period,  e.g.  due  to  loop  closure,  will  be  seen  there. 

Appendix  C  lists  the  short-period  approximation  transfer  functions  of  the  plants  actually 
used  in  the  CAD  package.  See  Section  4.4  for  more  detail  on  the  selection  of  these  plants. 

4.2  Design  Technique 

The  effective  plant  model  used  in  the  QFT  design  includes  not  only  aircraft  dynamics,  but 
also  actuator  and  sensor  dynamics.  Actuator  models  are  those  given  in  Eq.  (2.8)  through  Eq.  (2.10). 
Swift  did  not  account  for  sensor  dynamics  in  his  Lambda  model.  Sensor  dynamics  are  therefore 
included  in  the  state-space  produced  by  his  model. 

Several  methods  exist  for  the  design  of  a  digital  controller.  In  the  first  method,  design  can  be 
performed  directly  in  the  z-plane.  This  method  is  not  frequently  used  due  to  numerical  accuracy 
problems  which  arise  in  the  mapping  of  the  s-plane  to  the  z-plane  via  the  transformation  z  = 
where  T  is  the  sampling  period.  Another  method  [9]  is  to  transform  the  z-plane  plant  matrix 
into  the  to'-plane  via  the  transformation  z  =  Lacey’s  Lambda  design  was  performed  in 

the  w'-plane.  Problems  arose  with  the  gain  introduced  by  the  w'  method.  Therefore,  the  Pseudo- 
Continuous  Time  (PCT)  method  is  utilized  for  this  design,  [9]. 
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Implementation  of  the  PCX  method  requintrb  the  dominant  poles  and  zeros  of  the  plant  to  lie 
"sufficiently  close”  to  the  origin.  For  a  pole  located  at  atom  4-  sufficiently  close  is  defined 

as  [9]: 


0.1 

T 

(4.1) 

l^dom  1  ^ 

0.6114 

T 

(4.2) 

For  a  sampling  rate  of  bQHz,  the  sampling  period  T  isT  =  =  -^  =  0.02sec.  Therefore,  to  apply 

the  PCX  approach: 


(4.3) 

<  30.57 

(4.4) 

This  requirement  is  satisfied  by  the  full  equations  for  all  longitudinal  plants.  Note  that  actuator 
dynamics  do  not  lie  within  this  region.  The  PCX  approximation  is  valid  since  actuator  dynamics 
are  not  dominant. 

The  PCX  method  allows  the  zero-order  hold  effects  at  the  input  to  the  plant  to  be  modeled 
using  the  first-order  Pade  approximation.  The  zero-order  hold  is  modeled  exactly  in  the  s-plane 
as  [9]: 

1  -  e-*^ 

G.okis)  =  - -  (4.5) 

s 

The  first-order  Pade  approximation  models  G^oh  linearly  as: 
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This  appnndmation  is  valid  for 


where  Ug  is  the  plant  cutoff  frequency  and  u,  is  the  sampling  frequency  which  in  our  case,  corre¬ 
sponds  to  50Hz  and  is  w,  —  lOOir  =  314^.  Thus,  from  Eq.  (4.7)  Wj  =  31,416^.  Figure  4.2  shows 
the  plant  frequency  response  with  actuator  dynamics  included.  It  is  seen  that  the  cutoff  frequency 
is  approximately  w  =  15^.  Thus,  Eq.  (4.7)  is  satisfied. 


Figure  4.2  Longitudinal  Plant  and  Actuator  Dynamics 


Sampler  attenuation  of  f  must  also  be  included  with  the  PCT  model.  The  linear  nature 
of  the  PCT  approximation  allows  the  sampling  effects  to  be  included  with  the  zero-order  hold. 
Therefore,  the  PCT  approximation  Ct^(s)  is  written  as: 


G>i(^-)=^G,,,(s) 


(4.8) 
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Therefore  combining  Eq.(  4.6)  and  Eki.  (4.8)  with  a  sampling  time  of  T=0.02  seconds  yields: 


100 

s  +  100 


(4.9) 


The  PCT  implementation  of  the  Lambda  lon^tudinal  channel  is  shown  in  Figure  4.3. 


4.3  Loaded  and  Effective  Plants 

The  equations  defining  the  effect  of  control  surface  failures  on  Lambda  dynamics  are  derived 
in  Chapter  3.  These  equations  are  implemented  in  the  modified  Lambda  model  macro.  Once 
implemented,  the  altered  macro  generates  A  and  B  matrices  for  both  the  longitudinal  and  lateral 
channels.  A  Matlab  macro  then  converts  the  A  and  B  matrices  for  each  plant  into  a  data  file 
compatible  with  Mathematica.  This  is  required  because  the  MIMO  QFT  CAD  Package  operates 
from  a  Mathematica  platform.  Once  loaded  into  Mathematica  and  the  MIMO  CAD  Package, 
the  state-space  representations  are  transformed  to  transfer  functions  relating  output(s)  to  input(s). 
QFT  design  is  performed  with  the  resulting  transfer  functions. 
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Two  sets  of  A  and  B  matrices  are  generated  by  Swift’s  amended  macro.  One  set  defines  the 
bare  aircraft  dynamics,  Pi(s).  Matrices  for  the  disturbance  dynamics  Pd(^)  are  also  created.  As 
shown  in  Chapter  HI,  the  A  matrix  for  Pi{s)  and  Pd(«)  is  the  same  for  the  longitudinal  design. 
However,  B  matrices  differ  by  a  factor  of 

The  longitudinal  design  is  a  MISO  design.  For  the  tracking  case,  there  is  one  input  qemd  &Dd 
one  output  q.  Pi{s)  therefore  consists  of  a  1  x  1  matrix  of  transfer  functions.  For  the  disturbance 
rejection  case,  there  is  one  input  dtrim  and  one  output  q.  Similarly,  Pd{s)  is  also  a  1  x  1  matrix  of 
transfer  functions.  These  transfer  functions  are  listed  in  Appendix  C.  Henceforth,  both  Pi(s)  and 
Pd(s)  for  the  longitudinal  channel  are  denoted  as  scalar  quantities.  Figure  4.4  shows  the  frequency 
response  of  P{(s)  and  Pd(s)  for  all  plant  cases  used  in  the  QFT  design. 


Figure  4.4  Frequency  Response  of  Longitudinal  Pj(s)  and  Pd{s) 
PI  Response  Pd  Response 


radfsec  rad/isec 
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Controller  synthesis  with  the  QFT  technique  takes  into  account  not  only  aircraft  plant  dy¬ 
namics,  but  also  actuator  and  sensor  dynamics.  Zero-order  hold,  actuator  and  plant  dynamics  are 
combined  into  an  effective  plant  transfer  function  P,.  Figure  4.3  shows  the  relation  of  Pe  to  the  QFT 
MISO  model.  Figure  4.5  shows  the  frequency  response  of  the  P^  used  in  the  design  (short-period 
P,).  Appendix  C  lists  the  effective  short-period  longitudinal  plants. 


Figure  4.5  Frequency  Response  of  Lon^tudinal  P,  (short-period) 


4.4  Plant  Templates 

Plant  templates  are  generated  from  the  effective  plants  Pe(^)’  A  number  of  frequencies  must 
be  chosen  to  use  for  template  generation.  Templates  below  u?  =  1^  are  not  used.  Plant  data 
below  this  frequency  is  not  valid  since  the  short-period  approximation  is  used.  The  upper  and  lower 
tracking  bounds  for  a  m.p.  system  are  shown  in  Fig.  4.6  for  the  longitudinal  tracking  specification 
given  in  Table  1.1.  The  desired  tracking  response  bandwidth  Wh„  is  defined  as  the  frequency  at 
which  the  —12dB  line  intersects  the  upper  tracking  bound  [8].  Figure  4.6  shows  the  desired  tracking 
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respcHue  is  negligible  above  Wk^  =  20^.  Therefore,  templates  for  frequencies  above  w  =  20^ 
are  not  critical  in  the  design  for  tracking  and  disturbance  purposes.  High  frequencies  however,  still 
must  be  considered  for  stability.  One  high  frequency  template  at  w  =  50^  is  used  to  shape  the 
high  frequency  U-contour.  In  between  w  =  1  and  w  =  20,  template  frequencies  are  chosen  on  the 
basis  of  plant  variation.  Where  plants  vary  significantly,  frequencies  are  chosen  closer  together.  For 
the  longitudinal  design,  the  frequencies  w  =  1,2.5,3.25,4,6,8, 10,20  and  50  are  used. 

Figure  4.6  Longitudinal  Closed- Loop  Tracking  Bounds 


Templates  are  generated  by  determining  the  magnitude  and  phase  of  the  response  for  each 
effective  plant  at  a  given  frequency.  The  template  is  defined  by  plants  with  the  outermost  magnitude 
and  phase  values.  Lines  drawn  between  these  plants  define  the  template  boundaries.  Figure  4.7 
shows  tl>  iongitudinal  templates  generated  by  the  MIMO  QFT  CAD  Package. 

To  reduce  package  computations,  plsmts  which  do  not  define  a  boundary  at  any  of  the  chosen 
frequencies  are  discarded.  This  process  reduces  the  number  of  plants  to  15. 
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Several  flight  conditions  not  among  the  initial  32  are  then  included  in  templates  to  see  if  the 
template  boundaries  are  expanded.  One  iteration  of  this  process  is  described.  It  is  seen  that  little 
phase  variation  between  plants  exists  in  the  high  and  low  frequency  templates.  For  this  reason, 
an  intermediate  frequency  template,  w  =  6,  is  analyzed.  Figure  4.8  shows  that  plants  2,3,7,8,5,14 
and  11  lie  consecutively  on  the  boundary.  Segments  between  each  of  these  plants  comprise  the 
boundary.  Seven  boundary  segments  exist  for  this  template. 

Plants  5  and  8  comprise  a  segment.  The  contribution  of  plant  4  to  that  boundary  is  considered 
negligible.  The  only  flight  condition  which  varies  between  plants  5  and  8  is  speed.  Adding  another 
plant  with  the  same  flight  conditions,  but  with  speed  chosen  half-way  between  that  of  plants  5 
and  8,  yields  plant  16  which  expands  the  original  template  boundary.  Figure  4.9  shows  the  new 
template  with  plant  16  added.  Addition  of  this  plant  results  in  a  more  curvi-linear  boundary. 
Similar  template  expansion  occurs  for  several  other  frequency  templates. 
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Similar  analysis  for  the  other  boundary  segments  found  plants  which  lie  inside  the  existing 
boundaries.  Generally,  the  templates  assume  a  "banana”  shape.  The  lower  boundary  (segment 
11-14)  of  the  template  thus  creates  some  overdesign  in  the  system. 


The  16  plants  used  in  the  design  are  listed  in  Table  4.1.  The  corresponding  short-period 
transfer  functions  are  the  ones  listed  in  Appendix  C. 


Thble  4.1  Sixteen  Final  Longitudinal  Flight  Conditions 


Plant# 

eg 

(in) 

Speed 

(ftlsee) 

Qoo 

Weight 

(lbs) 

C 

1 

45.83 

118.23 

12.27 

215 

1 

2 

45.83 

185.79 

39.83 

215 

1 

3 

45.83 

185.79 

30.31 

215 

1 

4 

47.75 

118.23 

16.13 

190 

1 

5 

47.75 

118.23 

12.27 

215 

1 

6 

45.83 

185.79 

39.83 

190 

1 

7 

45.83 

185.79 

39.83 

215 

1 

8 

45.83 

185.79 

30.31 

215 

1 

9 

45.83 

118.23 

16.13 

190 

.75 

10 

45.83 

118.23 

12.27 

215 

.75 

11 

45.83 

185.79 

39.83 

215 

.75 

12 

45.83 

185.79 

30.31 

215 

.75 

13 

47.75 

118.23 

16.13 

190 

.75 

14 

47.75 

118.23 

12.27 

190 

.75 

15 

47.75 

185.79 

39.83 

190 

.75 

16 

47.75 

152.00 

20.29 

190 

1 

The  maximum  magnitude  differential  in  a  template  is  16.6  dB  at  w  =  10.  The  Tna.YitnnTn 
phase  difference  is  80  degrees  at  w  =  6.  Investigation  of  the  location  of  plants  within  each  template 
reveals  that  speed  and  altitude  cause  the  largest  plant  variation  with  a  constant  percentage  of 
elevator  damage.  Elevator  damage  causes  a  loss  of  magnitude  and  a  slight  increase  in  phase  lag. 

The  QFT  procedure  requires  the  selection  of  a  nominal  plant.  This  plant  is  the  one  actually 
used  during  the  controller  synthesis  (loop-shaping)  process.  Tracking  and  disturbance  bounds  on 
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the  Nichd’s  Chart  axe  created  from  the  nominal  plant.  The  nominal  plant  is  usually  chosen  as 
the  plant  with  the  largest  phase  lag  and  smallest  magnitude  for  the  most  frequencies.  Plant  14 
of  Table  4.1  is  chosen  as  the  nominal.  It  is  seen  that  Plant  14  is  at  the  lowest  speed,  highest 
altitude  and  maximum  fiulure.  It  is  expected  that  control  efforts  will  be  most  difficult  under  these 
conditions.  The  nominal  plant  transfer  function  is: 

M.  =  4.5870(3  +  2.1700) 

Seis)  (3+ 1.3559  ±i2.9567)  ^  ’  '' 

4.5  Tracking  and  Disturbance  Bounds 

Figures  4.10  and  4.11  shows  the  tracking  and  disturbance  bounds  generated  by  the  MIMO 
QFT  CAD  Package.  Tracking  bounds  are  comparable  to  those  generated  in  a  Lambda  longitudinal 
design  by  Wright  Laboratories.  Tracking  bounds  are  created  from  the  templates  shown  in  Fig.  4.7. 
Disturbance  bounds  are  created  based  on  a  disturbance  rejection  limit  of  0.2.  That  is,  the  maximum 
system  output  from  a  disturbance  should  not  exceed  20%  of  the  disturbance  input.  Because  no 
specifications  for  disturbance  rejection  are  given,  0.2  is  arbitrarily  chosen.  If  actual  responses  are 
unsatisfactory  or  if  a  0.2  bound  cannot  be  met,  the  bounds  can  be  changed  accordingly. 

A  low  system  bandwidth  is  beneficial  to  feedback  design.  Large  bandwidths  increase  control 
effort  as  well  as  allow  the  introduction  of  sensor  noise  and  structural  dynamics.  In  a  digital  design, 
low  system  bandwidth  is  especially  crucial,  as  a  high  bandwidth  requires  faster  sampling  rates.  It  is 
seen  that  satisfying  disturbance  bounds  requires  an  excessively  high  bandwidth.  The  advantage  of 
QFT  is  evident  here.  To  minimize  bandwidth,  engineering  judgement  is  used  to  selectively  violate 
bounds.  Bounds  at  frequencies  with  minimal  effect  on  system  response  aire  ignored.  More  effort  is 
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jdaced  on  meeting  the  low  frequency  bounds.  This  is  because  higher  frequencies  experience  greater 
attenuation  by  the  plant. 


Figure  4.10  Longitudinal  Tracking  Bounds 
TK 


4.6  Compensator  Gn{s) 

The  longitudinal  compensator’s  transfer  function  Gn(s)  is  obtained  by  shaping  the  nominal 
loop  on  the  Nichol’s  Chart  so  that  all  frequencies  lie  above  their  corresponding  bounds  and  outside 
their  corresponding  stability  contour. 

Elevator  failures  increase  plant  phase  lag  and  decrease  plant  magnitude.  Both  factors  tend 
to  increase  the  bandwidth  of  the  open-loop  plants.  This  causes  a  significant  amount  of  lead  and 
gain  to  be  placed  in  the  controller.  Poles  of  the  compensator  must  therefore  be  placed  far  from  the 
ori^n,  so  as  not  to  negate  the  benefits  of  the  zeros  which  are  added.  This  creates  aliasing  problems 
when  the  controller  is  implemented  digitally.  The  poles  of  (7ii(^)  must  therefore  be  ’’slow”  enough 


80  Gii(2)  at  hdHz  emulates  Gn(«).  A  rule  of  thumb  to  be  foUowed  while  loop-shaping  is  to  keep  the 
natural  frequency  of  the  poles  less  than  |  the  sampling  frequency  (a;„  <  At  SOifz,  therefore, 
poles  should  be  placed  at  s  >  —150  as  a  minimum.  This  is  as  non-dominant  as  the  poles  can  be 
made  when  including  poles  to  obtain  a  proper  order  compensator. 

Problems  arise  when  loop-shaping  with  compensator  poles  >  -150.  The  nominal  loop  trans¬ 
mission  tends  to  violate  the  U-contour  at  high  frequencies.  Anticipating  the  additional  lag  of  the 
compensator  poles,  as  the  nominal  loop  is  designed,  a  phase  margin  of  60”  is  targeted.  Then,  when 
compensator  poles  are  placed,  the  additional  lag  of  these  poles  returns  the  nominal  loop  closer  a 
45”  phase  margin. 

Even  with  the  phase  margin  considerations  while  loop-shaping,  U-contour  violations  occur  at 
high  frequencies.  It  is  decided  that  some  low  frequency  boundaries  are  to  be  violated  to  complete 


the  loop  shape.  Tracking  bounds  are  not  as  restrictive  and  they  can  ail  be  met  fairly  easily. 
Disturbance  bounds  are  therefore  violated. 

Initial  efforts  at  loop-shaping  result  in  a  compensator  which  satisfies  all  tracking  bounds  but 
violates  all  disturbance  bounds  except  w  =  1.  Simulations  revealed  excellent  tracking  response  and 
disturbance  rejection.  However,  analysis  of  sensor  noise  effects,  (see  Sec.  4.10),  revealed  significant 
amplification  of  sensor  noise  between  5  and  90^.  This  noise  amplification  should  be  attenuated 
to  2u:ceptable  levels  as  aircraft  structural  modes  lie  within  the  specified  bandwidth.  Therefore,  a 
new  Gii(s)  is  synthesized. 

The  requirement  for  an  acceptable  G(s)  is  it  must  have  as  low  as  possible  magnitude  in  the 
bandwidth  of  concern.  This  presents  a  problem  with  respect  to  stability  and  the  need  to  provide 
quite  a  bit  of  additional  phase  lead  and  gain  due  to  failures  in  this  region.  Thus  it  is  decided 
that  loop-shaping  should  be  performed  to  satisfy  only  the  tracking  bounds.  The  first  loop-shaping 
attempt  violated  disturbance  bounds  significantly,  yet  disturbance  rejection,  upon  simulation,  was 
quite  good.  Ignoring  disturbance  bounds  allows  the  compensator  gain  to  be  lowered  considerably, 
thereby  enhamcing  noise  rejection. 

Loop-shaping  results  in  the  nominal  open-loop  transmission  of  Fig.  4.12  and  the  following 
s-domain  compensator: 


13804.6(s  +  2)(s  +  3)(s  +  13.125  ±  j7.26) 
s{s  +  0.3)(5  -I-  130)(5  -I- 150)2 


(4.11) 
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Based  on  the  tracking  bounds  in  Fig.  4.12,  the  initial  attempts  at  prefilter  design  indicate  that 
the  tracking  bound  at  w  =  2^  is  not  satisfied.  To  meet  the  tracking  bound  at  w  =  2^  requires  a 
significant  increase  in  compensator  gain  if  compensator  order  is  not  increased.  It  is  decided  that  the 
increase  in  compensator  order  is  not  worth  the  improvement  in  tracking.  Therefore,  the  bound  at 
ti7  =  2  is  left  unsatisfied.  With  this  violation,  simulation  will  determine  if  the  design  is  acceptable. 

Figure  4.12  Lon^tudinal  Nominal  Loop  Transmission 


The  digital  flight  control  system  of  the  Lambda  requires  the  controller  to  be  implemented 
discretely.  Because  Gnis)  does  not  have  two  more  poles  than  zeros  and  does  involve  a  zero- 
order-hold,  a  direct  z- transform  cannot  be  performed  [9].  Instead  the  following  bilinear  Tustin 
transformation  is  used  to  transform  (7ii(s)  to  Gix(z): 


2(1 -z-i)]’ 

r(i-bz-i). 


(4.12) 
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Applying  the  Tustin  approximation  to  Eq.  (4.11)  results  in  the  following  discrete  controller: 


12.25156(z  +  1)(2  -  0.96078)(z  -  0.94175)(z  -  0.76070  ±  jO.11300) 
{z  -  l){z  -  0.99402)(2  +  0.2)2(2  +  0.13043) 


(4.13) 


The  discrete  controller  must  adequately  duplicate  the  continuous  compensator  for  the  given 
sampling  rate.  The  duplication  can  be  determined  by  analyzing  the  frequency  responses  of  both 
Gii{s)  and  Gii(z).  If  the  responses  are  the  same  up  to  the  sampling  rate  is  adequate.  For 
u,  =  314.16^,  the  controller  frequency  responses  must  match  up  to  105^.  Figure  4.13  shows  a 
good  correspondence  within  the  required  bandwidth  of  0  <  w  <  105^. 


Figure  4.13  G'ii(s)  and  Git(z)  Controller  Frequency  Responses 

•  aiiM  -.-.•aii(z) 


4- 7  Prefilter  Fii(s) 

The  prefilter  is  obtained  by  plotting  the  frequency  responses  of  the  closed-loop  system  on  the 
same  magnitude  Bode  plot  as  the  upper  and  lower  tracking  bounds.  The  prefUter  is  then  adjusted 
until  the  frequency  responses  for  all  plants  lie  within  these  tracking  bounds.  If  violations  of  the 
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tracking  bounds  occur  in  the  loop  shaping  process,  it  is  likely  that  the  upper  and/or  the  lower 
tracking  bounds  will  be  violated  by  the  closed-loop  response  for  one  or  more  of  the  plants  at  that 
same  frequency.  If  tracking  bound  violations  occur  at  significantly  small  magnitudes,  the  impacts 
are  negligible.  The  prefUter  is  chosen  as: 


0.010942(5-1- 1.9)(s -I- 100)2 
(5-|-2.2)(5-H6.3)(5-|-15) 


(4.14) 


The  non-dominant  zeros  at  s  =  —100  are  placed  to  prevent  jPn('Z)  from  having  two  zeros  lie  on 
the  unit  circle,  as  would  occur  if  -f\i(s)  is  left  as  a  first  over  third  order  transfer  function  and  the 
Tustin  approximation  applied  to  obtain  Fx^{z). 

Figure  4.14  shows  that  with  the  prefilter  of  Eq.  (4.14),  the  tracking  bounds  are  slightly 
violated  for  some  plants  at  approximately  2^.  This  violation  is  expected  due  to  the  violation  of 
the  w  —  2  tracking  bound.  Once  again,  the  gain  required  to  meet  the  tracking  specifications  at 
these  frequencies  is  not  considered  worthy  of  the  cost  in  noise  attenuation. 

Like  <Jii(s),  i^ii(5)  must  be  transformed  to  the  discrete  domain.  A  Tustin  transformation 
results  in: 

^  0.035698(z-  0.96271)(z)2 

^  (2 -0.95695)(z-0.88147)(«- 0.73913)  ^  ‘  ' 

Since  the  i^u(5)  and  Fii(^)  frequency  responses  in  Fig.  4.15  lie  essentially  on  top  of  one 
another  within  the  bandwidth  0  <  «;  <  105  then  the  prefilter  of  Eq.  (4.14)  is  valid  for  the  given 
sampling  rate. 


4-20 


Figure  4.14  Longitudinal  Closed  Loop  Response 


Figure  4.15  i^ix(«)  and  Fii(z)  Pre-Filter  Frequency  Responses 

-F11W 
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4.8  Validation  of  the  QFT  Design 


A  simulation  of  the  QFT  design  reveals  if  any  of  the  design  specifications  are  violated.  It 
has  already  been  shown  in  Fig.  4.14  that  the  frequency  domain  tracking  specifications  are  slightly 
violated  with  the  designed  controller  and  prefilter.  The  simulation  (Sec.  4.9)  will  determine  if  this 
violation  is  insignificant.  Wright  Laboratories  also  specified  a  7  =  45”  phase  margin  angle.  This 
corresponds  to  a  3dB  Mm  contour  on  the  NichoFs  Chart.  Figure  4.16  shows  that  the  7  requirement 
is  satisfied  for  all  plants. 

Figure  4.16  Longitudinal  Open-Loop  Transmissions  with  Gii(s)  Included 


Open-Loop  Phase  (deg) 


Analysis  of  the  closed-loop  disturbance  response  in  Fig.  4.17  indicates  that  the  arbitrarily 
chosen  0.2  disturbance  bound  is  violated  from  2  to  35—.  The  simulations  will  determine  if  this 
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violation  is  acceptable.  The  violation  is  most  likdy  to  be  unacceptable  if  the  input  disturbance 
signal  has  a  large  magnitude  at  one  of  the  violated  frequency  components. 

Figure  4.17  Longitudinal  Closed-Loop  Disturbance  Response 


rad^ 


4.9  Simulation 

Time  domain  simulations  verify  the  effectiveness  of  the  designed  contrcriler.  Simulations  are 
performed  in  Simulink  using  a  linear  model.  A  non-linear  actuator  element  is  included  in  the 
model.  The  non-linear  element  represents  the  deflection  limits  of  the  elevator.  SimtUink  allows 
both  discrete  and  continuous  elements  to  be  included  in  the  simulation.  Therefore,  prefllters  and 
controllers  are  described  in  the  model  with  z-domain  transfer  functions  with  sampling  time  equal 
to  0.02  seconds.  Actuators  and  aircraft  dynamics  are  described  with  s-domain  transfer  functions. 
Figure  4.18  is  the  block  diagram  representation  of  the  simulation  model. 
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Figure  4.18  Longitudinal  Channel  Simulink  Model 


A  step  input  of  q-ref  =  is  applied  to  the  system.  Figure  4.19  shows  the  aircraft  response 
with  the  designed  controller  and  prefilter  in  place.  The  response  is  seen  to  reach  the  desired  reference 
in  a  reasonable  amount  of  time.  Pitch  angle  increases  steadily  as  expected.  The  corresponding 
elevator  deflection  is  well  within  the  limits  of  the  actuator. 

The  controller  is  seen  to  attenuate  the  pitch  motion  caused  by  the  loss  of  elevator  and  also 
return  the  aircraft  to  trimmed,  level  flight.  The  elevator  is  seen  to  assume  a  steady-state  value. 
This  elevator  setting  is  a  result  of  the  trim  moment  lost  when  the  elevator  damage  occurs. 

A  larger  reference  input  of  10^  is  applied  to  the  system  to  analyze  possible  saturation  of 
the  actuators.  Figure  4.20  shows  the  elevator  saturating  as  early  as  1.25  seconds  for  one  flight 
condition.  The  flight  condition  which  saturates  (flight  condition  #10)  is  characterized  by  low  speed 
and  high  altitude  with  elevator  failure.  Minimal  control  authority  is  expected  here.  When  the 
elevator  saturates,  the  pitch  rate  is  not  m2dntsuned.  This  is  due  to  the  corresponding  loss  of  speed. 
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When  actually  performing  a  pitch  maneuver,  the  pilot  must  apply  more  throttle  to  maintain  pitch 
rate.  The  response  of  Fig.  4.20  is  considered  adequate  for  the  Lambda  aii craft. 

Figure  4.19  Longitudinal  IVacking  and  Disturbance  Responses 


There  exists  an  extent  of  damage  at  which  no  controller  can  adequately  maintain  mrcraft 
dynamic  response.  Obviously,  if  Q  =  0,  no  capability  exists  to  control  the  longitudinal  dynamics 
of  the  aircraft.  While  the  designed  controller  and  prefilter  of  Eq.  (4.13)  and  Eq.  (4.15)  generally 
meet  the  given  specifications  for  failures  up  to  Ce  =  0.75,  it  is  desired  that  the  controller  still  allow 
some  control  authority  for  more  drastic  surface  damage. 

When  flying  at  steady  cruise  conditions,  the  elevators  are  set  at  a  small  positive  deflection 
to  trim  the  aircraft.  As  surface  area  is  lost,  the  devators  deflect  more  to  generate  the  same  trim 
pitch  moment  on  the  tail.  There  exists  a  failure  level  when  the  damaged  elevator  saturates  before 
creating  the  necessary  trim  moment.  When  this  occurs,  the  aircraft  is  statically  unstable.  If  the 
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Longitadinal  Tracking  Response  10^  Reference 
Tnckn  AnBonit 


Oi  1  2  2i 


Oi  1  1i  2  2i  3 


elevator  does  trim  the  aircraft  but  is  saturated,  no  control  authority  is  left  in  the  elevator.  It  is 
therefore  desirable  to  maintain  some  available  elevator  deflection  for  dynamic  control  purposes. 

The  total  elevator  deflection  after  damage  is  equal  to  the  initial  trim  setting  (^tr<m)  plus  the 
additional  deflection  required  for  damaged  trim  (£').  To  maintain  5”  of  elevator  control  authority, 
the  following  relation  must  hold: 


Strin,  +  <  10“ 


(4.16) 


Simulations  assume  a  worst  case  initial  trim  setting  of  5”.  Therefore,  the  steady-state  elevator 
deflection  in  response  to  a  longitudinal  disturbance  should  not  exceed  5”  to  satisfy  £q  (4.16). 
Equation  (3.75)  defines  the  relation  between  Strim  and  tf'.  Substituting  Eq.  (3.75)  into  Eq.  (4.16) 


Assuming  the  worst  case  initial  trim  setting  of  5*,  solving  £q.  (4.17)  for  Ce  results  in: 

C.  >  0.5  (4.18) 

Therefore,  the  elevator  cannot  experience  more  than  a  50%  failure  and  still  guarantee  elevator 
control  authority.  Figure  4.21  shows  a  Lambda  response  for  the  16  flight  conditions  with  Ce  =  0.5. 


Figure  4.21  Longitudinal  Response  C«  =  0.5 
Tnckira  RMpoflM  OittffeMMRMpOM 


It  is  seen  that  the  additional  trim  deflection  is  equal  to  5**.  The  elevator  deflection  in  response 
to  a  1^  pitch  rate  reference  never  exceeds  5”.  Therefore  reasonable  longitudinal  control  is  possible 
with  only  50%  of  available  elevator.  The  rate  response  though  is  not  as  crisp  as  at  75%  available 
elevator.  This  is  because  QFT  design  was  based  on  25%  elevator  failure.  Larger  magnitude  of  pitch 
references  saturate  the  elevator.  These  simulations  are  based  on  am  initial  worst  caise  trim  setting 
of  5**.  This  is  a  conservative  estimate. 
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4.10  Other  Concerns 


The  Lambda  has  longitudinal  structural  modes  at  28^.  Excitation  of  the  elevator  by  a 
reference  input  or  by  sensor  noise  should  avoid  exciting  these  modes.  Block  diagram  models  in 
Simulink  are  used  to  determine  the  frequency  response  of  the  elevator  with  respect  to  both  of  the 
cited  inputs. 

Figure  4.22  shows  the  equivalent  block  diagram  used  for  analysis  and  Fig.  4.23  shows  the 
resulting  frequency  response  of  a  reference  input.  It  is  seen  that  the  system  is  greatly  attenuated 
above  10^.  This  greatly  eliminates  the  probability  of  pilot  induced  oscillations  (p.i.o.). 

Figure  4.22  Block  Diagram  of  Elevator  Response  to  Reference  Input 


Aicnft 


Figure  4.24  shows  the  equivalent  block  diagram  used  for  analysis  and  Fig.  4.25  shows  frequency 
response  of  an  elevator  to  sensor  noise.  Noise  levels  are  approximately  7  to  8dB  lower  than  in  the 
initial  compensator  design.  While  sensor  noise  is  not  amplified,  the  possibility  exists  that  if  sensor 
noise  is  large  enough,  structural  mode  excitation  may  become  significant. 
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Figure  4.23  Frequency  Response  of  Elevator  to  a  Reference  Input 


w 


Figure  4.24  Block  Diagram  of  Elevator  Response  to  Sensor  Noise 


Son  V  lef 
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Figure  4.25  Frequency  Response  of  Elevator  to  Sensor  Noise 


The  effect  that  pitch  rate  sensor  noise  has  on  elevator  control  can  be  determined  through 
simulation.  Wright  Laboratories  has  developed  white  Gaussian  noise  models  for  all  sensors  on 
board  the  Lambda  [10].  The  modeled  pitch  rate  sensor  noise  is  added  to  the  Simulink  model  at 
the  sensor  feedback  loop.  Figure  4.26  shows  that  pitch  rate  noise  has  a  negligible  effect  on  the 
aircraft’s  pitch  channel  response. 
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V.  Lateral  MIMO  Design 


This  chapter  presents  the  MIMO  lateral  controller  design.  The  design  of  the  lateral  channel 
differs  slightly  from  the  design  of  the  MISO  longitudinal  channel.  The  differences  are  discussed 
here.  Individual  design  steps  are  outlined  in  Chapter  1. 


5.1  Plant  Generation 

5.1.1  Flight  Conditions.  Similar  to  the  longitudinal  channel,  flight  conditions  used  to 
generate  lateral  channel  plants  are  selected  to  adequately  represent  the  entire  flight  envelope.  The 
same  flight  variables  in  the  longitudinal  channel  determine  the  envelope  for  the  lateral  channel; 
center  of  gravity,  airspeed,  altitude  and  weight.  The  same  range  of  these  variables  used  in  the 
longitudinal  channel  is  used  to  generate  lateral  plants.  Elevator  failure  is  still  considered  in  the 
lateral  channel  due  to  the  cross-coupling  of  failures.  Failures  to  the  lateral  control  surfaces,  mlerons 
and  rudders,  are  now  considered.  This  leaves  seven  different  variables  to  be  considered  for  the 
lateral  channel.  Varying  all  seven  between  their  maximum  and  minimum  values  results  in  128 
permutations  of  lateral  plants.  Processing  this  many  permutations  during  the  design  process  would 
be  computationally  intensive.  The  number  of  plants  is  therefore  reduced. 

Because  the  lateral  channel  has  no  gravity  force  acting  on  it,  (see  Eqs.  (3.27)  through  (3.29)), 
weight  has  no  effect  on  lateral  dynamics.  Inspection  of  the  magnitude  and  phase  response  of 
the  lateral  transfer  functions  reveals  that  the  lateral  response  is  negligibly  affected  by  the  center  of 
gravity  location.  These  two  variables  are  therefore  set  to  the  median  value  of  their  expected  ranges: 


c.g.  =  45.83”  -|- 


47.74”  -  45.83” 
2 


46.79” 


(5.1) 
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215/63  -  190/6s 

weight  =  190  + - - - = 


(5.2) 


The  nmnber  of  lateral  plants  is  reduced  further  by  assuming  a  constant  elevator  failure  of  50% 
(C«  =  0.5).  A  const  tralue  is  assumed  because  the  magnitude  of  elevator  failure  has  no  effect  on 
lateral  disturbance  dynamics.  It  only  effects  the  magnitude  of  the  disturbance  input.  Therefore,  if 
disturbances  from  50%  elevator  failure  are  rejected,  failures  of  less  magnitude  are  also  rejected. 

The  four  parameters  that  are  varied  for  lateral  plant  generation  are; 

•  Airspeed 

•  Altitude 

•  Aileron  Failure 

•  Rudder  Failure 

which  result  in  sixteen  lateral  flight  conditions.  These  fliight  conditions  are  listed  in  Appendix  B. 
Like  the  longitudinal  channel,  plants  are  eliminated  during  template  generation.  This  is  discussed 
further  in  Section  5.4. 


5.1.2  Plant  Transfer  Function  Matrix.  Lateral  control  system  design  can  take  place 
around  several  different  lateral  states.  Because  coordinated  flight  is  desired,  roll  rate  and  sideslip 
angle  are  the  states  utilized  for  design.  No  approximations  are  used  to  generate  these  transfer 
functions. 


The  resulting  transfer  function  matrix  utilizing  bare  aircraft  dynamics  is: 


p(»)  p(») 

«.(•)  <r(.) 

PM  iisl 
«.(.)  SrM 


(5.3) 
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5.2  Design  Technique 


The  PCT  method  is  used  for  the  lateral  design.  Dominant  poles  of  aircraft  plant  transfer 
functions  all  lie  within  the  region  specified  by  £q.  (4.3)  and  Eq.  (4.4).  Figure  5.1  shows  the 
frequency  response  of  lateral  plant  and  actuator  dynamics.  Cutoff  frequencies  satisfy  the  limit  of 
30^  defined  by  Eq.  (4.7).  The  PCT  approximation  is  therefore  valid  for  the  lateral  channel. 


Figure  5.1  Lateral  Plant  and  Actuator  Dynamics 
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The  PCT  implementation  of  the  lateral  channel  is  shown  in  Fig.  5.2.  Gx{s)  is  given  by 
Eq.  (4.9).  Actuator  dynamics  are  given  in  Biqs.  (2.9)  and  (2.10).  W  =  {«>,;}  is  a  2x2  weighting 
matrix  which  is  designed  as  an  aileron- rudder  interconnect. 
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Figure  5.2  Lateral  PCT  Implementation 


5.S  Loaded,  Effective  and  MISO  Equivalent  Plants 

The  roll  and  yaw  channels  are  designed  sequentially.  For  presentation  purposes,  the  designs 
are  described  simultaneously  here. 


5.S.1  Loaded  Plants.  Like  the  longitudinal  channel,  lateral  plants  are  generated  with 
the  modified  Matlab  macro  in  Appendix  A.  Three  different  matrices  are  generated  for  the  lateral 
channel.  One  A  matrix  for  each  flight  condition  is  generated.  This  matrix  is  the  same  for  both 
the  tracking  and  disturbance  dynamics.  Two  B  matrices  are  created,  one  for  Pj(s)  and  one  for 
Pd(s).  From  the  appropriate  A  and  B  matrices,  the  MIMO  CAD  Package  creates  the  2x2  tracking 
dynamics  matrix  P/(5)  and  the  2x  1  disturbance  dynamics  matrix,  Pd(s).  Transfer  functions  related 
to  the  /3  output  contain  non-minimum  phase  (n.m.p.)  elements.  Pi(s)  is  shown  in  Eq.  (5.3).  Pd(s) 
is  shown  as: 


P^(^) 


p<(«) 

.(» 


•(» 


(5.4) 
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The  transfer  functions  which  comprise  Pj(s)  and  for  ^  flight  conditions  used  in  the 

lateral  design  are  listed  in  Appendix  C.  The  frequency  responses  of  Pi(s)  and  Prf(s)  are  shown  in 
Fig.  5.3  and  Fig.  5.4  respectively. 


Figure  5.3  Lateral  Pi{s)  Frequency  Responses 


Figure  5.4  Lateral  P4(«)  Frequency  Responses 


5.S.2  Effective  Plants.  Bare  aircraft  dynamics  are  combined  with  actuator  dynamics, 
PCT  dynamics  and  a  weighting  matrix  to  form  the  effective  plants  Pe(s).  These  plants  are  used 
in  the  QFT  loop-shaping  process. 

The  2x2  weighting  matrix  W  is  used  as  an  aileron-rudder  interconnect.  When  Lambda  rolls 
in  response  to  an  mleron  deflection,  sideslip  is  induced.  A  non-zero  element  in  W2,i  causes  the 
rudder  to  deflect  whenever  the  aileron  deflects.  The  induced  rudder  deflection  counters  the  natural 
tendency  of  the  aircraft  to  generate  sideslip  when  responding  to  a  roll  command.  The  reduction 
of  the  output  due  to  a  roll  command  by  use  of  a  weighting  matrix  minimizes  the  control  efforts 
required  by  the  QFT  controller.  Likewise,  a  non-zero  Wi^2  element  causes  the  aileron  to  deflect 
when  a  command  is  sent  to  the  rudder.  Because  the  rudder  is  not  intended  to  be  used  often  in  a 
coordinated  aircraft  and  because  the  roll  induced  by  a  rudder  command  is  negligible,  u;i,2  is  set  to 
0.  Diagonal  terms  are  set  equal  to  1  since  no  ampliflcation  or  attenuation  of  command  signals  is 
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desired.  To  keep  the  system  order  to  a  minimum,  W2,i  is  considered  as  only  a  gain  element.  No 
additional  dynamics  are  introduced. 


The  weighting  matrix  element  W2,i  is  chosen  through  analysis  of  the  aircraft  open-loop 
time  response.  The  open-loop  response  includes  actuator  dynamics.  The  value  of  W2,i  which 
best  minimizes  for  a  short  time  span  is  1.0.  Figure  5.5  compares  the  open-loop  response  of 


W  = 


Wi.i  «?i.j 

W2,i  tV2,7 


1  0 
0  1 


with  W  = 


1  0 
1  1 


Figure  5.5  Weighting  Matrix  Effect  on  Open-Loop  Response 
8oadUra->w(2,1H>.0  DalidLkw->w(2,1H0 
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It  is  seen  that  when  W2,i  =  1.0  the  median  excursion  of  /?  is  centered  about  0  =  0.  The  rudder 
induces  a  positive  yaw  rate  to  maintain  0  as  close  to  0  as  possible.  An  adverse  effect  of  a  non-zero 
element  in  W2,i  is  an  increase  in  the  peak  the  0  response  due  to  a  disturbance  input.  A  large  pitch 
rate  is  induced  by  the  modeled  lateral  disturbance.  The  aileron  command  to  counter  that  large 
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rate  interconnects  to  the  rudder  causing  the  rudder  to  over  compensate  yaw  motion.  However,  the 
postive  results  of  cross-coupling  rejection  outwdgh  the  negative  effects  of  disturbance  rejection. 

frequency  responses  of  the  effective  plants  with  W21  =  1.0  are  shown  in  Fig.  5.6.  Appendix 
C  lists  the  effective  lateral  plants. 


Figure  5.6  Lateral  Pe(s)  Frequency  Responses 


5.S.S  MISO  Equivalent  Loops.  The  MIMO  QFT  method  for  an  mxm  plant  requires 
the  formation  of  m*  equivalent  MISO  loops  for  design  as  shown  in  Fig.  5.7  [5].  Compensators 
are  designed  in  the  same  manner  as  a  simple  1x1  QFT  design.  Each  MISO  loop  has  a  tracking 
input  and  a  disturbance  input.  The  tracking  input  is  the  external  reference  command  applied  to 
the  plant.  The  disturbance  input  represents  the  cross-coupling  from  other  channels.  The  diagonal 
compensator  and  prefilter  matrices  are  given  by: 


gn  0  0 

G(s)  =  0  fl22  0  (5.5) 

0  0  P33 

■  /u  0  0 

F(s)  =  0  /22  0 

0  0  /as 

Therefore  the  off-diagonal  loops  only  have  disturbance  inputs. 


The  equivalent  loop  transfer  functions  qu  are  formed  by  taking  the  inverse  of  the  plant  matrix 
Pe(s)  and  then  inverting  each  element  of  P7*(s). 


Pn  P23 
Ph  P*33 


and: 

1  1 

Q(s)=  ""  ”” 

1  1 
.  Pm  Pm 

The  elements  of  Q(s)  are  the  qa  transfer  functions  which  appear  in  the  equivalent  loops  for  each 
plant  case. 
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Two  methods  exist  for  forming  the  equivalent  loops.  The  first  method,  Method  One,  treats 
the  loops  as  entirely  separate  QFT  designs.  Loop  transmissions  Li  from  previously  completed  loops 
are  not  considered  in  cross-coupling  disturbance  rejection.  The  second  method,  Method  Two,  is 
an  improvement  on  the  first  in  that  it  takes  into  account  loop  transmissions  of  previously  designed 
loops.  Previously  completed  loops  alter  the  disturbance  input.  When  starting  the  design  process. 
Method  One  is  used  for  loop  2.  For  loop  3,  Method  Two  is  used  which  requires  the  generation  of 
a  new  expression  for  933  [11].  This  new  expression  is  denoted  as  933^  and  is  given  by: 


?33«  —  933 


1  +  La 

1  +  Lj  —  7j3 


(5.9) 


where: 


P37P73 

P33P32 


(5.10) 


A  general  rule  of  thumb  exists  for  choosing  the  first  loop  to  be  designed.  Normally,  the  loop 
with  the  smallest  bandwidth  is  designed  first  which  for  this  design  is  the  yaw  channel.  However, 
because  turn  coordination  is  desired,  the  roll  channel  is  designed  first.  By  designing  the  roll  channel 
first,  roll  loop  transmissions  are  modeled  in  the  yaw  channel.  This  allows  loop-shaping  to  precisely 
eliminate  cross-coupling  of  roll  commands  into  as  desired  for  turn  coordination. 

The  MIMO  CAD  Padu^e  forms  the  Q(s)  matrices  from  the  Pe(s)  matrices  and  in  turn  forms 
933e.  Frequency  response  data  of  933(3)  and  qa3t{s)  are  shown  in  Fig.  5.8. 
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Figure  5.7  m’  Ek^uivalent  MISO  Loops 


Figure  5.8  Frequency  Responses  of  922(«)  (left)  and  qssei*)  (right) 


5.4  Plant  Templates 


A  set  of  templates  is  generated  for  both  lateral  channels  from  the  frequency  response  of  932(5) 
933«(^)-  Plant  template  frequencies  are  chosen  as  described  in  Sec.  4.4.  For  the  lateral  channel, 
frequencies  of  w  =  0.5, 1,2, 3.5, 5, 8, 10,20  and  100  are  chosen. 


As  in  the  longitudinal  channel,  plants  which  do  not  define  a  template  boundary  in  either 
loop  are  eliminated.  Two  plants  are  also  added  to  expand  the  boundaries.  The  design  is  therefore 
accomplished  with  thirteen  plants.  Table  5.1  lists  the  thirteen  flight  conditions  used  in  the  lateral 
design.  Plant  templates  generated  with  the  thirteen  flight  conditions  are  shown  in  Fig.  5.9  and 
Fig.  5.10. 


Table  5.1  Thirteen  Final  Lateral  Flight  Conditions 


# 

eg 

(in) 

Speed 

ift/sec) 

(»»//*’) 

Weight 

ilbs) 

c. 

Cr 

c. 

1 

46.79 

118.23 

16.13 

202.5 

0.5 

1 

0.5 

2 

46.79 

118.23 

12.27 

202.5 

0.5 

0.5 

1 

3 

46.79 

185.79 

39.83 

202.5 

0.5 

1 

1 

4 

46.79 

185.79 

39.83 

202.5 

0.5 

1 

0.5 

5 

46.79 

185.79 

30.31 

202.5 

0.5 

0.5 

1 

6 

46.79 

185.79 

30.31 

202.5 

0.5 

0.5 

0.5 

7 

46.79 

118.23 

16.13 

202.5 

0.5 

1 

1 

8 

46.79 

118.23 

16.13 

202.5 

0.5 

1 

0.5 

9 

46.79 

118.23 

12.27 

202.5 

0.5 

0.5 

1 

10 

46.79 

185.79 

39.83 

202.5 

0.5 

1 

1 

11 

46.79 

185.79 

30.31 

202.5 

0.5 

0.5 

1 

12 

46.79 

152.00 

26.66 

202.5 

0.5 

1 

0.5 

13 

46.79 

152.00 

20.29 

202.5 

0.5 

1 

1 

Examination  of  plant  location  within  the  templates  of  each  channel  reveals  that  failures  cause 
a  drop  in  gain  and  a  slight  increase  in  phase  lag.  In  addition  to  failures,  gain  variations  are  caused 
by  velocity  changes.  Gain  increases  with  velocity.  Velocity  and  altitude  variation  cause  the  greatest 
variation  in  phase.  The  maximum  variation  of  template  parameters  are: 
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Figure  5.9  Lateral  Plaat  Templates  -  RoU  Chauinel 
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Figure  5.10  Lateral  Plant  Templates  -  Yaw  Channel 
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•  Loopl  -  17  dB  @  «;  =  100 

24"  ®w  =  5 

•  Loop2  -  21  dB  @  u;  =  100 

106*  @  u;  =  3.5 


Plant  ^6  is  chosen  as  the  nominal  plant  for  both  channels.  The  nominal  plant  transfer 
functions  are: 


5.5  Loop-Shaping 

Tracking,  cross-coupling  and  disturbance  rejection  bounds  are  generated  with  the  MIMO 
QFT  CAD  Package.  Like  the  longitudinal  channel,  not  all  bounds  can  be  met  and  still  allow 
a  reasonable  system  bandwidth.  Loop-shaping  is  therefore  performed  with  emphasis  on  meeting 
tracking  bounds.  High  frequency  bounds  are  readily  violated  to  mmntain  a  low  system  bandwidth. 

Figures  5.11  and  5.12  show  the  nominal  open-loop  transmissions  for  channel  1  and  2  respec¬ 
tively  with  Gii^a)  =  1.  Tracking  bounds  are  displayed.  It  is  seen  that  bounds  for  the  yaw  channel 
are  much  less  restrictive.  Initially,  bounds  at  =  0.5^  are  not  included  in  the  design  process 
which  results  in  significant  tracking  errors  in  the  yaw  channel  over  the  time  span  of  interest.  While 
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tracking  in  the  yaw  channel  is  not  of  critical  interest,  it  is  still  desirable  to  have  some  degree  of 
control  in  that  channel  with  0  feedback  employed.  Addition  of  a  bound  at  u;  =  0.5^  and  adjust¬ 
ment  of  the  loop  shape  greatly  improves  tracking.  The  improved  tracking  also  improves  disturbance 
rejection  because  of  the  increased  gain.  Tracking  performance  can  be  improved  further,  but  it  is 
desired  to  keep  a  relatively  low  gain  and  bandwidth  for  the  yaw  channel  due  to  the  noise  associated 
with  feedback. 

5.6  Compensators  Gii{s) 

Lambda  lateral  channels  have  lower  bandwidths  than  the  longitudinal  channel.  This  greatly 
facilitated  loop-shaping  as  poles  of  the  lateral  compensators  are  not  as  difficult  to  place.  Loop¬ 
shaping  results  in  the  following  compensators: 


and: 


^22(5)  = 


6200(s  +  0.1)(s  -f-  3)(s  +  9±  >15.59) 
s\s  +  125)(s  -I-  130)(s  -I- 150) 


(5.13) 


-5250(5  +  0.75)(s  +  1.6)(s  1.5  ±  >2) 

s{s  +  0.1)(s  -I-  15)(3  -I-  45)(s  +  50) 


(5.14) 


The  physical  implications  of  the  URV  response  require  that:  (1)  the  yaw  channel  loop  trans¬ 
mission  have  a  -180°  shift  in  phase,  thus  6^33(5)  must  have  a  negative  gain;  and  (2)  the  yaw  rate  r 
and  sideslip  /9  have  opposite  signs  in  the  gain  term  of  their  respective  transfer  functions  which  also 
requires  that  the  feedback  signals  have  their  sign  inverted. 
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Figures  5.13  and  5.14  show  the  final  loop-shape  tor  the  roll  and  yaw  channel  respectively. 
Note  the  tracking  violations  at  low  frequencies  for  the  roll  channel.  These  violations  are  assumed  to 
have  negligible  effect  on  the  tracking  response.  Time  domain  simulations  validate  this  assumption. 


Application  of  the  Tustin  approximation  to  £q.  (5.13)  and  £q.  (5.14)  results  in  the  following 
x-domain  controllers: 


5.99046(z  +  1)(«  -  0.99800)(x:  -  0.94175)(s  -  0.79808  ±  j0.25717) 
~  {z-  l)’(z  +  0.2)(2:  +  0.13043)(2:  +  0.11111) 


and: 


,,  -22.12114(z+ l)(z-0.98511)(z- 0.96850)(«- 0.96968  ±0.038811) 

~  (z  -  l)(z  -  0.99800)(z  -  0.73913)(z  -  0.37931)(z  -  0.33333) 
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Figure  5.14  Yaw  Channel  Final  Loop-Shape 


Comparison  of  the  s-domain  and  a-domain  frequency  responses  of  each  controller  show  good 
correlation  over  the  region  of  interest.  See  Fig.  5.15. 


I 


5.7  Prefilters  Fii{9) 


The  following  prefilters  are  formed  for  each  channel; 


0.006(s-t- 100)* 
~  (s-|-6)(s-|-10) 


(5.17) 


and: 


0.0008(s+l)(s-|-100)* 

(s  -I-  2)(s  +  1.4  ±  j 1.42829) 


(5.18) 


Similar  to  the  longitudinal  prefilter,  nondominant  zeros  are  added  to  the  lateral  prefilters. 
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Figure  5.15  Lateral  ControUer  frequency  Responses 


-022(1}  •.•.-Q22(;9  -033(1)  -.-.-033(4 


Figure  5.16  shows  the  closed-loop  Requency  response  for  each  tracking  response  and  for  each 
cross-coupling  response.  It  is  seen  that  roll  channel  tracking  bounds  are  violated  above  w  =  10^. 
This  corresponds  to  the  violation  of  tracking  bounds  on  the  Nichol’s  Chart.  Violations  do  not  occur 
until  the  closed-loop  response  is  significantly  attenuated.  The  roll  channel  closed-loop  frequency 
response  is  therefore  considered  adequate. 

Tracking  bound  violations  also  occur  in  the  yaw  channel.  Initial  design  attempts  result  in  a 
controller  and  prefilter  which  meet  the  specified  yaw  tracking  bounds.  However,  the  gain  of  this 
controller  is  large  enough  to  make  the  implementation  of  it  impractical.  The  high  gain  controller 
amplifies  sensor  noise  and  causes  significant  chatter  in  the  rudders.  A  new  controller  is  therefore 
designed  which  somewhat  violates  tracking  bounds,  but  does  not  excessively  amplify  sensor  noise 
(see  Sec.  5.10).  The  tracking  bounds  could  be  better  satisfied  if  a  higher  order  prefilter  is  used. 
Because  of  sampling  rates  and  the  extent  to  which  the  yaw  channel  is  used  for  control,  a  higher  order 
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piefilter  is  not  deured  either.  Simulations  determine  the  adequacy  of  the  designed  yaw  channel 
controller. 


Cross-coupling  rejection  is  seen  to  be  quite  good.  Violations  occur  in  the  roll  response  due 
to  a  rudder  input.  The  violations  are  not  corrected  because  the  rudder  is  not  a  primary  command 
surface  and  the  violations  do  not  result  in  inadequate  time  domain  responses. 

The  z-domain  representation  of  the  lateral  compensators  derived  via  the  Tustin  transforma¬ 
tion  is: 

F33(z)  =  0.020583; - — T  (5.19) 

^  ^  (2  -  0.88679)(z  -  0.81818)  '  ^ 

olid: 

^33(^)  =  0.0030811t - n  Q<in7fi\7  4.  0007777 

'  (2 -0.96078)(z-l- 0.97200  ±j0.027777  ^  ’ 

Figure  5.17  shows  the  2-domain  representations  are  valid  over  the  region  of  interest. 

5.8  Validation  of  the  QFT  Design 

It  is  shown  that  the  closed-loop  tracking  responses  are  adequate  in  meeting  the  given  specifica¬ 
tions.  Open-loop  transmissions  shown  in  Fig.  5.18  and  Fig.  5.19  show  that  the  7  =  45°  specification 
is  satisfied  in  both  channels.  Figure  5.20  shows  that  the  chosen  disturbance  responses  are  also  sat- 


OpeivLoop  Ptwe  (ctod) 


S.9  Simtdation 


Time  domain  simulations  verify  the  effectiveness  of  the  designed  lateral  controllers.  A  Simulink 
model  with  non-linear  actuator  and  rate  gyro  elements  is  used  for  the  simulation.  As  in  the  lon¬ 
gitudinal  channel,  prefilters  and  controllers  are  implemented  in  the  z-domain  while  aircraft  and 
actuator  dynamics  are  in  the  s-domain.  Figure  5.21  shows  the  Simulink  block  diagram  used  for 
the  simulation  of  the  lateral  channel. 

Figure  5.21  Lateral  Channel  Simulink  Model 


p_gyro 


A  unit  step  input  is  applied  to  each  of  the  control  inputs  at  different  times.  Figure  5.22  shows 
the  lateral  response  to  a  1^  roll  input.  It  is  seen  that  the  roll  rate  assumes  the  reference  value  very 
quickly.  Also,  the  cross-coupling  is  limited  very  effectively.  The  peak  value  of  /3  is  only  0.075®  for 
this  input.  Most  plants  remain  even  closer  to  zero.  It  is  seen  that  elevator  and  rudder  deflections 
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are  minimal.  Figure  5.23  shows  the  response  to  a  reference  command  of  1‘.  Note  that  the  ^ 
response  is  not  as  "tight”  as  the  p  response.  This  is  because  of  less  stringent  specifications  on  the 
response  since  yaw  is  not  the  primary  control  channel.  The  cross-coupling  into  the  roll  channel  is 
seen  to  be  minimized  by  the  controller.  The  maximum  value  of  the  roll  angle  is  only  0.08®.  Control 
surface  deflections  are  well  within  actuator  limits. 

Figure  5.22  Roll  Channel  Tracking  Response 

o. - - - 


3 
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Figure  5.24  shows  the  lateral  response  to  an  elevator  deflection.  The  elevator  deflection 
assumes  that  one  entire  elevator  is  lost.  The  remaining  elevator  surface  is  deflected  to  its  maximum 
value  of  15®.  The  roll  rate  is  seen  to  be  effectively  attenuated  by  the  controller,  limiting  the  roll 
angle  displacement  to  a  maximiun  value  of  0.9®.  The  aileron  is  seen  to  assume  a  non-zero  steady- 
state  value.  This  deflection  is  required  to  counter  the  roll  moment  induced  by  the  damaged  elevator. 
The  maximum  steady-state  deflection  of  -3®  still  allows  room  for  the  aileron  deflection  to  maintain 
control  authority  in  the  roU  axis.  The  maximum  value  of  the  sideslip  angle  is  0.8®  in  the  disturbance 


response.  This  value  can  be  reduced  by  lowering  the  magnitude  of  w^^i.  However,  the  negative 
effects  of  W3,i  on  disturbance  are  far  outweighed  by  the  gains  in  cross-coupling  rejection. 


Actuator  and  gyro  saturation  effects  are  investigated  only  for  the  primary  input  channel  (roll). 
A  reference  step  input  of  40^  is  applied  to  the  system.  It  is  seen  that  the  roll  rate  starts  to  run 
away.  This  is  due  to  saturation  in  the  roll  gyro.  Therefore,  for  a  50%  aileron  failure,  the  gyros  will 
saturate  first.  See  Fig.  5.25. 


5.10  Other  Concerns 

As  in  the  longitudinal  channel,  pilot  oscillations  (p.i.o.)  and  sensor  noise  are  concerns.  The 
frequency  response  of  three  transfer  functions  are  analyzed:  roll  reference  to  aileron,  roll  gyro  noise 
to  aileron  and  P  sensor  noise  to  aileron.  The  transfer  function  is  not  analyzed  since  yaw  control 
inputs  are  avoided  with  this  design. 
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Figure  5.24  Lateral  Channel  Disturbance  Response 
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Figtires  5.26  through  5.29  show  the  equivalent  block  diagrams  of  the  aileron  output/reference 
input  signal  and  the  aileron  output/sensor  noise  input  signal  configurations  and  the  corresponding 
frequency  responses.  It  is  seen  that  high  frequency  pilot  oscillations  and  roll  gyro  noise  are  well 
attenuated. 

Figure  5.26  Block  Diagram  of  Aileron  Output /Reference  Input 


■mfi 


The  effect  of  noise  is  a  concern  in  the  yaw  channel  especially  when  0  feedback  is  to  be  used. 
The  P  sensor  is  a  noisy  sensor  as  it  measures  unsteady  airflow  dynamics.  Excessive  amplification  of 
sensor  noise  can  cause  "chatter”  in  the  rudders.  Figure  5.30  shows  the  equivalent  block  diagram  of 
the  rudder  response  to  sensor  noise  input  and  Fig.  5.31  shows  the  frequency  response  of  the  aileron 
to  sensor  noise  transfer  function.  Amplification  of  the  noise  appears  to  be  significant.  However, 
the  frequency  response  with  the  chosen  Ggg,{s)  is  significantly  lower  than  the  initial  yaw  channel 
design.  Amplification  of  the  noise  above  5^  is  due  to  high  gain  in  the  G33  controller  element  (see 
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Figure  5.27  Frequency  Response  of  Aileron  to  Reference 


Figure  5.28  Block  Diagram  of  Aileron  Output/Roll  Sensor  Noise  Input 
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Figure  5.29  Frequency  Response  of  Aileron  to  Sensor  Noise 


Fig.  5.15).  When  surface  failures  occur,  the  plant  transfer  functions  experience  a  loss  of  gain  and 
additional  phase  lag.  To  counter  these  effects,  gain  and  phase  lead  are  returned  to  the  open-loop 
transmission  with  the  compensator.  This  additional  gain  and  phase  lead,  while  restoring  desired 
open-loop  qualities,  also  amplifies  sensor  noise.  The  redesigned  controller,  while  not  perfectly 
meeting  yaw  channel  specifications,  maintains  sensor  noise  amplification  at  reasonable  leveb. 

Time  domain  simulations  are  performed  using  white  Gaussian  noise  models  of  the  Lambda 
sensors.  Figure  5.32  and  Fig.  5.33  show  that  efforts  to  reduce  noise  effects  on  the  control  surface 
outputs  have  been  successful.  Figure  5.32  shows  the  time  domain  response  of  the  lateral  channel  to 
roll  rate  gyro  noise.  Roll  gyro  noise  is  seen  to  have  a  negligible  effect  on  the  aircraft  states.  Both 
the  aileron  and  rudder  experience  a  ±0.1”  in  response  to  the  noise.  The  rudder  chatters  due  to  the 
aileron-rudder  interconnect.  Figure  5.33  shows  the  time  domain  response  of  the  lateral  channel  to 
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Figure  5.30  Block  Diagram  of  Rudder  Output //9  Sensor  Noise  Input 
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sensor  noise.  Again,  effects  on  the  lateral  states  are  negligible.  The  rudder  experiences  a  ±0.2* 
chatter.  The  aileron  chatter  is  extremely  small  since  there  is  no  rudder-mleron  interconnect. 


Figure  5.32  Response  of  the  Lambda  to  Roll  Sensor  Noise 


If  farther  redaction  of  noise  effects  in  either  channel  is  desired,  the  controller  can  be  redesigned 
for  a  smaller  degree  of  damage  to  the  control  surfaces.  This  eases  phase  and  gain  requirements  on 
the  controllers,  thus  minimizing  noise  amplification.  Another  approach  can  be  to  include  a  sixth- 
order  Butterworth  filter  in  the  feedback  loops  to  filter  sensor  noise  [1].  This  needs  to  be  included  at 
the  start  of  the  QFT  design  process  to  ensure  robust  behavior  of  the  resulting  controller.  Another 
option  is  to  design  the  QFT  controller  around  roll  rate  p  and  yaw  rate  r.  The  yaw  gyro  is  not 
expected  to  be  as  noisy  as  the  /?  sensor.  However,  coordination  efforts  will  not  be  as  direct  as 
designing  with  jd  feedback. 
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VL  Autopilot  Design 


This  chapter  presents  an  autopilot  design  for  the  Lambda.  The  autopilot  is  designed  around 
the  QFT  p,  0  controller  synthesized  in  the  previous  chapters.  The  designed  QFT  controller  performs 
the  functions  of  minimizing  sideslip  and  stabilizing  flight  over  a  wide  range  of  parameters.  While 
roll  rate  is  an  excellent  state  to  design  a  QFT  controller  around,  it  is  desired  to  allow  the  pilot 
to  command  a  desired  bank  angle.  The  commanded  bank  angle  determines  how  tight  a  turn  the 
aircraft  will  make. 

The  autopilot  is  a  very  simple  design.  It  simply  integrates  the  output  of  the  roll  gyro  to  yield 
the  bank  angle  0.  This  design  is  therefore  limited  by  the  saturation  rates  of  the  roU  gyro,  ^40^. 
The  bank  angle  calculated  by  the  int^ator  is  compared  to  a  reference  bank  angle  input  by  the 
pilot  to  determine  the  reference  error.  The  pilot  also  commands  a  pitch  rate  to  initiate  the  banked 
turn.  The  pitch  rate  and  sideslip  angle  commands  are  input  directly  to  the  QFT  controller. 

Figure  6.1  shows  the  Simulink  model  used  to  simulate  the  autopilot.  Both  the  longitudinal 
and  lateral  channels  are  simulated  simultaneously  to  see  the  effect  of  damaged  elevator  deflections 
on  the  lateral  channel.  The  model  includes  the  discretized  controllers  and  preiilters,  continuous 
actuator  and  plant  models,  and  all  nonlinear  elements.  The  thirteen  lateral  flight  conditions  are 
used  since  they  contain  the  most  drastic  elevator  damage.  These  conditions  best  analyze  the 
rejection  of  elevator  damage  coupling  into  the  lateral  channel. 

The  first  simulation  shows  the  effect  of  damaged  elevator  deflections  on  the  lateral  channel. 
The  elevator  gradually  increases  its  deflection.  It  is  seen  that  for  all  plants  the  QFT  controller  and 
autopilot  eliminate  motion  in  the  lateral  channel  in  response  to  damaged  elevator  deflections  while 
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Figure  6.1  Autopilot  Simuliok  Model 


maintaining  the  commanded  pitch  rate.  Lateral  moments  generated  by  the  elevator  are  countered 
by  the  aileron  and  rudder. 

The  second  simulation,  Fig.  6.3,  shows  the  aircraft  response  to  a  45*  roll  command.  The 
aileron-rudder  interconnect  in  the  QFT  controller  maintains  adequate  turn  coordination.  The  peak 
value  of  P  is  about  3*.  However,  the  turn  coordinator  quickly  returns  to  0*.  The  maximum 
excursions  occur  for  the  simulation  performed  at  the  slowest  speed,  maximum  altitude  and 
maximum  damage  conditions,  as  expected.  This  is  precisely  the  flight  condition  with  the  least 
control  authority  .  However,  the  exhibited  control  is  still  adequate.  The  desired  roll  angle  is 
achieved  rapidly  and  precisely  for  all  flight  and  damage  conditions. 

Figure  6.4  shows  the  aircraft  response  to  a  roll  doublet  maneuver.  A  10*  doublet  of  2  sec¬ 
ond  duration  (each  pulse)  is  commanded.  The  roll  channel  is  seen  to  respond  in  a  smooth  and 
coordinated  manner.  Note  the  rudder  response  as  it  maintains  coordination. 
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The  next  simolation,  Fig.  6.5,  applies  a  10^  pitch  rate  reference  to  the  system.  All  flight 
conditions  have  a  lateral  channel  response  since  the  flight  conditions  used  always  assume  a  50% 
elevator  damage.  Tom  coordination  is  maintained  with  minimal  lateral  control  surface  deflection. 
As  shown  in  Chapter  IV,  the  commanded  pitch  rate  is  not  maintained  when  the  elevator  saturates. 
Saturation  occurs  more  quickly  than  in  previous  simulations.  During  actual  flight  the  pilot  can 
minimize  the  loss  of  pitch  rate  by  appl3dng  throttle  commands  to  maintain  airspeed. 

A  10^  pitch  doublet  of  2  second  duration  (each  pulse)  is  now  commanded.  Figure  6.6  shows 
the  response.  Pitch  rate  tracks  the  reference  adequately  until  the  elevator  saturates.  The  aileron 
and  rudder  are  seen  to  maintain  lateral  coordination. 

Although  rudder  control  is  not  intended  to  be  the  primary  control  method,  rudder  inputs 
are  simulated.  Figure  6.7  shows  the  aircraft  response  to  a  5*  sideslip  command.  This  amounts  to 
placing  the  aircraft  into  a  flat  turn.  The  desired  sideslip  angle  is  maintained.  Notice  the  aileron 
response  as  it  holds  the  wings  at  a  level  condition. 

Finally,  a  5”  sideslip  doublet  is  commanded  as  shown  in  Fig.  6.8.  The  aircraft  is  seen  to 
adequately  track  the  commanded  reference  signal. 
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VII.  Summary,  Conclusions  and  Recommendations 

This  chapter  presents  a  summary  of  the  thesis,  draws  conclusions,  and  recommends  future 
research. 

7.1  Summary 

The  Quantitative  Feedback  Theory  design  technique  is  applied  to  the  Lambda  URV  to  synthe¬ 
size  a  robust  flight  control  system.  This  design  expands  on  previous  Lambda  flight  control  system 
designs  in  that  it  considers  actuator  failures  as  a  source  of  structured  plant  uncertainty.  Failures 
considered  are  losses  of  surface  area  inflicted  on  the  elevators,  mlerons  and  rudders.  Mathematical 
models  of  the  failure’s  effects  on  the  dynamics  are  derived  and  added  to  an  existing  macro  which 
generates  state-space  models  of  the  Lambda  for  given  flight  conditions.  Failures  manifest  them¬ 
selves  as  changes  in  stability  derivatives  and  as  external  disturbances  in  both  the  longitudinal  and 
lateral  channels.  The  longitudinal  channel  design  thus  requires  a  MISO  controls  design  with  one 
feedback  variable.  The  lateral  channel  entails  a  MIMO  design  consisting  of  two  feedback  states. 
Finally,  an  outer-loop  autopilot  is  designed  around  the  QFT  controller. 

The  longitudinal  MISO  QFT  design  begins  with  the  selection  of  flight  conditions  which  ade¬ 
quately  represent  the  Lambda  flight  envelope.  Sixteen  flight  conditions  are  chosen.  The  short-period 
approximation  of  longitudinal  aircraft  dynamics  is  utilized.  The  PCT  design  technique  is  employed 
for  a  digital  sampling  rate  of  50Hz. 

In  the  pitch  channel,  a  fifth-order  controller  and  third  order  prefilter  are  synthesized  which 
employ  the  feedback  state  (pitch  rate)  to  adequately  meet  frequency  domain  stability  and  tracking 
specifications  in  the  face  of  varying  flight  conditions  and  elevator  failure.  The  controllers  are 
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transformed  to  the  z-plane  via  a  Tostin  transformation.  The  digital  design  is  simulated  with 
Simulink.  Analyses  of  the  possibility  of  P.I.O.S  and  sensor  noise  effects  are  then  performed. 

The  lateral  MIMO  QFT  design  proceeds  in  a  similar  manner.  Thirteen  flight  conditions 
are  chosen.  A  PCT  design  is  performed  around  two  feedback  states,  roll-rate  and  sideslip  angle. 
Sideslip  angle  is  chosen  as  a  feedback  state  so  the  resultant  QFT  controller,  which  in  a  natural  way 
renders  a  degree  of  decoupling,  coordinates  aircraft  flight.  To  aid  coordination,  an  aileron-rudder 
interconnect  is  employed  via  a  control  weighting  matrix.  An  equivalent  m  x  m  MISO  loop  transfer 
function  matrix  Q  is  developed  which  results  in  a  fifth-order  roll  channel  controller  and  a  second- 
order  prefllter.  Modified  yaw  channel  Q’s  are  formed  using  QFT’s  Method  Two.  In  addition, 
a  fifth-order  yaw  controller  and  a  third-order  prefilter  which  adequately  meet  frequency  domain 
specifications  are  synthesized.  Digital  simulations  of  roll  and  yaw  channel  disturbance  responses 
are  performed  with  Simulink.  Analyses  of  lateral  channel  P.LO.s  and  sensor  noise  effects  are 
performed. 

An  outer-loop  autopilot  is  formed  around  the  QFT  controUer.  The  autopilot  integrates  roll 
rate  to  obtain  roll  angle.  Roll  angle  is  then  fed  back  to  the  pilot  reference  command.  The  remote 
pilot  can  then  perform  a  banked  turn  by  simply  commanding  a  reference  bank  angle  and  a  pitch 
rate.  No  rudder  control  is  necessary  since  the  QFT  controller  induces  the  turn  by  eliminating 
sideslip.  The  autopilot  is  simulated  with  Simulink.  The  Simulink  model  includes  cross-coupling 
of  the  lon^tudinal  channel  into  the  lateral  chsmnd  as  derived  for  elevator  failures.  Seven  different 
commanded  inputs  are  simulated. 
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l.t  Conclusions 


Several  condosions  are  reached  about  this  specific  design  and  the  QFT  method  in  general: 

•  A  robust  digital  controUer  and  prefilter  are  designed  for  the  longitudinal  channel  which  meet 

frequency  domain  specifications  for  a  loss  of  elevator  surface  area  up  to  25%.  Reasonable 
aircraft  control  can  be  maintained  for  up  to  50%  elevator  loss  but  response  specifications 
are  no  longer  met.  Greater  damage  leveb  cannot  be  adequately  accommodated  due  to  the 
inherently  large  bandwidth  of  the  lon^tudinal  channel.  The  controller  quickly  returns  the 
aircraft  to  a  trimmed  flight  condition  after  surface  area  Ic  The  designed  controller  minimizes 

effects  of  pitch  rate  sensor  noise  on  the  time  responses. 

•  Robust  digital  controllers  and  prefilters  are  designed  for  the  lateral  channel.  The  resulting 
controllers  and  prefilters  adequately  meet  specifications  for  50%  surface  area  loses  in  both 
the  aileron  and  rudder.  The  controUers  reject  cross-coupling  of  the  longitudinal  channel 
into  the  lateral  channel  in  the  face  of  elevator  damage.  The  controllers  and  aileron-rudder 
interconnect  successfully  limit  sideslip  induced  by  an  aircraft  roll.  IVacking  specifications  are 
slightly  violated  because  of  controller  gain  limitations  caused  by  sensor  noise. 

•  An  autopilot  is  designed  which  enables  coordinated  flight  for  the  Lambda  URV  under  serious 
damage  conditions.  The  QFT  controller  and  the  autopilot  allow  the  Lambda  to  be  flown  in 
a  ’’fee^nc  the-floor”  manner.  That  is,  banked  turns  require  only  a  roU  angle  command  and  a 
pitch  rate  command.  Autopilot  performance  is  limited  by  the  aerodynamic  properties  of  the 
aircraft  itself  and  by  limits  on  the  gyro  feedback  rates. 

•  Aerodynamic  surface  failures  cause  significant  gain  and  phase  losses  in  the  aircraft  response. 
Therefore,  for  responses  to  meet  specifications,  the  controller  must  contain  significant  gain 
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aod  phase  lead.  Noise  and  sampling  effects  limit  the  ability  of  the  controllers  to  have  these 
characteristics.  Noise  limits  the  amount  of  gain  which  may  be  placed  in  the  controllers. 
This  is  especially  a  concern  when  using  a  state  such  as  0  whose  sensor  is  characteristically 
noisy.  Sampling  rate  limits  the  bandwidth  of  the  controller  which  may  be  used.  More  drastic 
failures  can  be  accommodated  using  a  QFT  controller  if  sensor  noise  is  decreased  and/or  a 
faster  sampling  rate  is  used. 

s  Control  surface  failures  can  be  successfully  included  as  structured  parameter  uncertainties  in 
a  QFT  flight  control  design.  QFT  robustly  compensates  for  limited  control  surface  failures 
to  maintmn  a  desired  aircraft  response.  The  robust  control  is  accomplished  with  a  fixed,  not 
a  scheduled,  contndler,  thus  maintaining  controller  simplicity.  The  utility  of  QFT  is  evident 
throughout  the  design.  Its  effective  use  of  the  frequency  domain  along  with  the  associated 
MIMO  QFT  CAD  Package  allows  quick  insight  into  the  flight  control  problem  when  design 
trade-offs  need  to  be  made. 

•  The  AFIT  MIMO  QFT  CAD  Package  is  invaluable  in  the  completion  of  this  and  any  QFT 
design.  Without  its  ease  of  use,  many  aspyects  of  this  design  problem  could  not  have  been 
investigated  as  effectively  or  as  eflSciently  as  they  were. 

7.S  Recommendations 

The  following  recommendations  are  made  based  upon  the  results  of  this  thesis: 

•  Future  study  should  investigate  the  improvement  of  the  aileron-rudder  interconnect  by  using 
dynamic  terms  in  the  weighting  matrix.  Better  rejection  of  sideslip  in  response  to  an  aileron 
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input  cnn  be  achieved  with  a  dynamic  matrix.  Improved  sideslip  rejection  will  in  turn  improve 
coordinated  flight  efforts. 

•  Greater  control  in  the  lonptudinal  channel  can  be  achieved  if  flaperons  are  used  for  dynamic 
contrcd.  Now  that  a  faulure  model  is  generated,  a  weighting  matrix  can  be  employed  to 
share  control  efforts  between  the  elevator /flaperon  and  aileron/flaperon  surfaces.  Differential 
surface  control  can  also  be  investigated  for  damaged  flight  control. 

•  Investigate  coordinating  damage  control  efforts  between  a  QFT  flight  control  system  and  a 
multiple  model  adaptive  estimator. 

•  Apply  QFT  to  design  a  flight  control  system  for  a  high  performance  Air  Force  aircraft  with 
control  surface  failures. 

•  Modify  MIMO  QFT  CAD  Package  to 
ticular,  options  allowing  the  placement 


allow  greater  freedom  when  generating  plots.  In  pair- 
of  labels  and  custom  text  on  plots  should  be  included. 
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Appendix  A.  Modified  Model  Generation  Macro 


The  following  Matlab  executable  file  was  created  by  Lt  Gerald  Swift  and  modified  by  Capt 
Steve  Rasmussen.  It  was  further  modified  to  account  for  changes  in  stability  derivatives  due  to 
control  surface  damage. 


This  Mcro  calculates  the  longitudinal  and  lateral- 
directional  state  space  equations  of  notion  for  the 
Lanbda  URV  given  leg,  Ul,  ql.  W  and  theta. 

Input :  X*[Xcg,Ul ,ql .H ,theta,ztae,ztar ,ztaa] 

share,  Xcg  «  inches 

Ul  •  feet/second 
ql  *  pounds/sq  foot 
H  •  pounds 
theta  *  degrees 

ztae  >  elevator  danage  (1  to  0;  !■  no  danage) 
ztar  •  rudder  danage  (1  to  0) 
ztaa  •  aileron  danage  (1  to  0) 


Outputs : 

Clong*  EClfu ,  CNa ,  CNad ,  CHq ,  CLu ,  CLa ,  CLad ,  CLq ,  CDa ,  Clhi , 

CLde,CDde,CHde,CLdf ,CDdf ,CMdf] ’ 

Clat> [Clb , Clp , Clr , Clda , Cldr , Cnb , Cnp , Cnr , Cnda , Cndr , 

Cyb,Cyp,Cyr,Cyda,Cydr] ’ 

Dlong  *  corresponding  vector  of  dinansional  derivatives  for  Clong 
Dlat  ■  corresponding  vector  of  dinansional  derivatives  for  Clat 
Along  ■  A  natrix  for  longitudinal  equations  of  notion 
Alat  «  A  natrix  for  lateral-directional  equations  of  notion 
Blong  ■  B  natrix  for  longitudinal  equations  of  notion 
Blat  ■  B  natrix  for  lateral-directional  equations  of  notion 
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t 

zcgvzd.l); 
a>z(1.2); 
q>z(1.3): 

■-z(1.4)/32.174: 
zt«««x(i,6) : 
zt«r"x(l,7}; 
ztaa*x(l,8): 

St»5.17; 

S«<tl .  1 : 
iia*.9S; 
k^.42: 

At«2.8645: 

AX-4.S5; 

dxcg* (46 . 405-xcg) / 12 ; 
lt-(101.99-xcg)/12; 

■■21 . 1 : 

€■1.51; 
b-14.07: 

Ixx«32.502: 

Iyj*26.666: 

IzzMO.939: 

Ixz*^.6922: 

CLl«x(l,4)/(q*s): 

CDl-0.027; 

CHl-O.O; 

Al«Ixz/Ixx: 

Bl«Ixz/Izz; 

ztaar  ■  ((ztaa-l)*1.4  ♦  (ztar-l)*1.0  +  24.6)/24.6;  t 


ztazpzB  ■  ((ztar-l)*1.0  *  3.5}/3.S;  X 
ztmapza  ■  ((zta«-l)*1.4  ■*'  5.17)/5.17;  X 
naa  ■  (l-zta«pza)«nu*St*At*lt*(l-k) ;  X 
daa  ■  S«*A«*dxcg  +  zta«pra*nu*St*At*lt*(l-k);  X 
daaaa  ■  naa/daa;  X 
nza  ■  (l>ztaepn)*nu*St*At*(l-k) ;  X 
dza  ■  Sv*Av  ztaapn*ntt*St*At*(l-k) ;  X 
dzaza  ■  nza/dza;  X 
X 


X  Tha  folloving  paraaatars  raprasant  changas  to 
X  Svifts’  iBodal  basad  on  20  nov  92  flight  taat  -  sjr 
X 


X  addad,  ask 
X  addad.  ask 
X  addadt  ask 
X  addad,  ask 
X  addad,  ask 
X  addad,  ask 
X  addad,  ask 
X  addad,  ask 
X  addad,  ask 
X  addad,  ask 
X  addad,  ask 


addad, 

addad, 

addad, 

addad, 

addad, 

addad, 

addad, 

addad, 

addad. 


ask 

ask 

ask 

ask 

ask 

ask 

ask 

ask 

msk 


CMAV-0.7S; 

CIIADV-0.7; 


CHqV«0.6; 

CIIDEV--0.5; 
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Cm?-0.7S; 

CDV-0.9: 

ciiiv-0.9: 

CIDRV—0.6; 

X 

X  «ind  nav  pariMtara  -  sjr 

X 

Clla<M}; 

Clto«CmV*(0.326643*z(l.l)-16.5189):  X  changad,  sjr 

Clfad-CIUDV*(0.238235*z(l.l}-7.8977):  X  chaogad,  sjr 

Cllq«ClfQV*(0.628420*x(l,l)<4S.0582);  X  changad,  ajr 

au«-0.0001719*x(l ,  1)40 .01518 : 

CLa«5.e20: 

CLad«-0.03838«x(l . 1)43.83508: 

CLq—0 .65202*x(l .  1)437 . 8884 ; 

COaaO; 

CDtt«0; 

CLda-0.2908: 

CDda^: 

Cflda<OfDEV«(0.013464*x(l.l)-1.48039):  X  chaogad.  sjr 
CLdf-1.419: 

CDdf *0.08489: 

Clldf*0.070932*x(l ,  l)-3.63838: 
Clongl*CCIfa.CHa,CNad,Cllq,CLu,CLa,CLad,CLq,CDa,CDvCI : 
Cloog2«CCLda.CDda.CHda,CLdf  .CDdl.CHdf] ; 

Cloiig*rcloiigl  Cloxig2]’: 

X 

Clb— 0.01451: 

Clp— .6538: 

Clr-0. 08763; 


Clda*CIiDAVaO . 2608 ;  X  chaogad.  uk 

Cldr-0 . 000213*x( 1 . 1) -0 . 00783 ; 

CDb«OBV*(-0.00038*x(l, 1)40. 07834);  X  chaogad.  sjr 

Cop*-0. 03601: 

Cor*CniV*(0.003087*x(l.l)-0. 31072);  X  chaogad.  sjr 

Cnda*-0. 01368; 

Cndr*CIDRV4(0.001633*x( 1,1) -0.15208);  X  updatad  4/22/91.  chaogad,  ajr 
Cyb*- 0.4372; 

Cyi>*-0.001600; 

Cyr—0 .00424*x(l ,  1)40 .46047 ; 

Cyda*0; 

Cydr-0.2865; 


Clat*[Clb, Clp, Clr.Clda.Cldr, Cob, Cop, Cor, Coda, Cndr, Cyb, Cyp.Cyr.Cyda.Cydr]  ' 
X 


X 
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X 

X 

X 

X 

X 

X 

X 


Calculation  of  Diaansional  DarlTativas 
Dlong«Otu,Ma,Mad,Nq.Zu,Za,Zad,Zq.Xa,Xu,Zde,X(le,Mde,Zdf  ,Xdf  ,N(lf] ' 
Dlat •  [Lb ,  Lp ,  Lr .  Lda ,  Ldr ,  Mb ,  Bp ,  Ir ,  Mda ,  Bdr ,  Yb ,  Yp ,  Yr ,  Yda ,  Ydr]  ' 


Hu»q*8*c*(CRu*2*Cl!l)/(Iyy*u) ; 

Ma"q*B*c<t‘(CMa-(CMa*dmama))/Iyy;  X  changed,  mak 

Nadaq*s*c*c*CMad/(2*l77*u) ; 

Nq*q*8*c*c*ztaepneCNq/(2*lyy*u) ;  X  changed,  msk 

Zu«-q*8*(CLu+2*CLl)/(m*u) ; 

Za*-q*8*((CLa-»C01)-'((CLa'»'CDl)edzaza))/n;  X  changed,  msk 
Zada-q*8*c*CLad/(2*m*u} ; 

Zq>-qe8ec*ztaeprmeCLq/(2*B*u) ;  X  changed,  msk 

Xu»-q*8*(CDu+2*CDl)/(m*u) ; 

Xaa-q*8* (CDa-CLl) /m : 

Zdea-q*8*ztae*CLde/m;  X  changed,  msk 

Zde«-q*8eCDde/m ; 

Mdeaq*8*c*ztae*CMde/Iy7:  X  changed,  msk 

Zdf  «-’qe8*CLdf /m : 

Xdf>-qe8*C0df/m: 

Ndf ■qKsoceCMdf / lyy : 

Dlong«Clfu,Ma,Mad,Nq,Zu,Za,Zad,Zq,Xa,Xu,Zde,Xde,Nde,Zdf  ,Xdf  ,Ndf]  ’ ; 
X 


Lb"qe8*beztarprmeClb/Izz ; 
Lpaq*s*beb*ztaar*Clp/(2*Ixz*u) ; 
Lraqe8*b*b*ztazpzm*Clr/(2*Ixz*u) ; 
Ldaaq*8*beztaa*Clda/Izz ; 
Ldr>q*8*b*ztar*Cldr/Ixz ; 
Bbaq*s*b*ztarpzB*Cnb/Izz ; 
Bpaq*s*beb*ztarpxiiieCnp/(2*Izz*u) ; 
Br*qe8eb*beztarpziB*Cnr/  (2*Izz*u) ; 
Bda>q*8abeztaa*Cnda/Izz ; 
Bdr>q*s*b«ztar*Cndx/Izz ; 
Yb*q*s*ztarpziB*Cyb/m ; 
Ypaqe8*b*ztarpxiB*Cyp/(2*m*u) ; 
Yr*q*8*beztarprm*Cyr/(2*m*u) ; 
Yda«q*8*zt aaeCyda/m ; 
Ydr*q*8*ztar*Cydr/m ; 


Dlat> [Lb , Lp , Lr , Lda , Ldr , Bb , Bp , Mr , Bda , Bdr , Yb 


changed,  msk 
changed,  msk 
changed,  lask 
changed,  msk 
changed,  msk 
changed,  msk 
changed,  msk 
changed,  msk 
changed,  msk 
changed,  msk 
changed,  msk 
changed,  msk 
changed,  msk 
changed,  msk 
,Yp,Yr,Yda,Ydr] ' ; 


X 

X  Calciilation  of  the  Longitudinal  A  and  B  matrices 
X 

theta-x ( 1 , 5) /57 , 29578 ; 
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Along«[Zu,Xa,0,-32. 174*co8(th«ta) ; 

Zu/(u-Zad) ,2a/(u-Zad) , (u+Zq)/(u-Zad) ,-32.174*Bin(thata) ; 

Nu+Nad^Zu/ (u-Zad)  ,Ma'*'Mad*Za/(u-Zad)  ,Mq+Mad*(u+Zq)/(u-Zad)  .... 
-Mad«32. 174*8in(thata) ; 

0,0,1,03; 

Blong*CZda;  X  flaps  deleted,  nsk 

Zde/(u-Zad} ; 

Mde+Mad*Zde/ (u-Zad) ; 

03; 

X 

X  Calculation  of  the  Lateral  Directional  k  and  B  matrices 
X 

BAel-Bl*!!; 

Alat»[Yb/u,Yp/u,32.174*coa(theta)/u,Yr/u-l;  X  psi  state  deleted,  msk 

(Lb*Al*Hb) /BA , (Lp+Al*Mp) /BA ,0 , (Lr+Al*Kr) /BA ; 

0. 1,0,0; 

(Mb+Bl*Lb)/BA, (Hp+BleLp)/BA,0, (■r+Bl*Lr)/BA3  ; 

Blat-CYda/u,Ydr/u;  X  psi  state  deleted,  msk 

(Ldaf Al*Hda) /BA , (Ldr+AleRdr) /BA ; 

0,0; 

(Rda-»^Bl*Lda)/BA,  (Hdr-fBl«Ldr)/BA3 ; 

X 

X  Generate  B  matrices  for  the  lateral  and  longitudinal  disturbances 
X 

blngdist  •  ((l>ztae)/ztae)*Blong; 

X 

Ldauf  ■  q*8*b*Clda/Izz; 

blatdist  ■  [0;((l>ztae}*Ldauf)/(6<fBA);0;(Bl*(l-ztae)eLdauf)/(6*BA)] ; 
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Appendix  B.  Flight  Conditions 


This  appendix  displays  the  flight  conditions  prior  to  template  generation  for  both  the  longi¬ 
tudinal  and  lateral  channels. 


B.l  Longitudinal  Conditions 


Table  B.l  Initial  Lon^tudinal  Flight  Conditions 


# 

eg 

(in) 

Speed 

(ft/sec) 

(16V/<*) 

Weight 

(lb*) 

c. 

1“ 

45.83 

118.23 

16.13 

190 

1 

2 

45.83 

118.23 

16.13 

215 

1 

3 

45.83 

118.23 

12.27 

190 

1 

4 

45.83 

118.23 

12.27 

215 

1 

5 

45.83 

185.79 

39.83 

190 

1 

6 

45.83 

185.79 

39.83 

215 

1 

7 

45.83 

185.79 

30.31 

190 

1 

8 

45.83 

185.79 

30.31 

215 

1 

9 

47.75 

118.23 

16.13 

190 

1 

10 

47.75 

118.23 

16.13 

215 

1 

11 

47.75 

118.23 

12.27 

190 

1 

12 

47.75 

118.23 

12.27 

215 

1 

13 

47.75 

185.79 

39.83 

190 

1 

14 

47.75 

185.79 

39.83 

215 

1 

IS 

47.75 

185.79 

30.31 

190 

1 

16 

47.75 

185.79 

30.31 

215 

1 

17 

45.83 

118.23 

16.13 

190 

0.75 

18 

45.83 

118.23 

16.13 

215 

0.75 

19 

45.83 

118.23 

12.27 

190 

0.75 

20 

45.83 

118.23 

12.27 

215 

0.75 

21 

45.83 

185.79 

39.83 

190 

0.75 

22 

45.83 

185.79 

39.83 

215 

0.75 

23 

45.83 

185.79 

30.31 

190 

0.75 

24 

45.83 

185.79 

30.31 

215 

0.75 

25 

47.75 

118.23 

16.13 

190 

0.75 

26 

47.75 

118.23 

16.13 

215 

0.75 

27 

47.75 

118.23 

12.27 

190 

0.75 

28 

47.75 

118.23 

12.27 

215 

0.75 

29 

47.75 

185.79 

39.83 

190 

0.75 

30 

47.75 

185.79 

39.83 

215 

0.75 

31 

47.75 

185.79 

30.31 

190 

0.75 

32 

47.75 

185.79 

30.31 

215 

0.7.5 

Appendix  C.  Longitudinal  Plant  Transfer  Functions 

This  appendix  lists  the  effective  plant  transfer  functions  and  the  disturbance  transfer 
functions  for  the  longitudinal  channel. 

C.l  Longitudinal  Peis)  Short-Period  Transfer  Functions 


Plant  1: 

qis)  85958.1000(5+2.1700) 

S,(s)  ~  is  +  1.3139  ±  j4.0857)(5+  10.8140  ±  j7.7387)(5  +  100) 

Plant  2: 

qis)  278729.4455(5  +  4.4800) 

^e(s)  ”  (s  +  2.7124  ±  i7.3369)(5  +  10.8140  ±  j7.7387)is  +  100) 

Plant  3: 

q(5)  212263.8319(5  +  3.4100) 

^e(5)  "  (5  +  2.0647  ±  j6.4204)(5+  10.8140  ±  ;7.7387)(5  +  100) 

Plant  4: 

qis)  109347.4836(5  +  3.0814) 

«e(5)  “  (5  +  1.8263  ±i3.5743)(5+  10.8140  +  j7.7387)(5+  100) 

Plant  5; 

qis)  83337.7877(5  +  2.0725) 

rfe(3)  ~  (5  + 1.2652  ±j3.1567)(5+ 10.8140  ±j7.7387)(5+ 100) 

Plant  6: 

9(5)  270021.2398(5  +  4.8422) 

«e(«)  ”  (s  +  2.8699  +  i5.6167)(5+  10.8140  +  j7.7387)(5  +  100) 

Plant  7: 

9(5)  270185.4425(5  +  4.2796) 

6eis)  ~  is  +  2.6118  +  j5.6715)(5  +  10.8140  +  j7.7387)(5  +  100) 

Plant  8: 

9(5)  205793.2696(5+3.2568) 

tfe(«)  ”  (s  +  1.9881  + j4.9606)(5+  10.8140  +  j7.7387)(5  +  100) 

Plant  9: 

9(5)  84612.1689(5+3.1857) 

de(«)  “  (5+ 1.8419 +  j4.4824)(5+ 10.8140 +  77.7387)(5+ 100) 
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Plant  10: 


g(3) _ 64468.5757(3  +  2.1433) _ 

6,(s)  ~  (3  +  1.2778  ±j3.9504)(«+  10.8140  ±  j7.7387)(3  +  100) 

Plant  11: 

9(3)  209047.0876(3+4.4250) 

«e(a)  “  (s  +  2.6377  ±  j7.0929)(3  +  10.8140  ±  j7.7387)(3  +  100) 

Plant  12: 

q{s)  159197.8739(3  +  3.3681) 

Se{s)  ~  (3  + 2.0079  ±i6.2078)(3+  10.8140  + j7.7387)(3+  100) 

Plant  13: 

9(3)  82010.6127(3  +  3.0415) 

Seis)  ~  (3  +  1.7809  ±j3.3700)(3+  10.8140  ±  j7.73S7)(3  +  100) 

Plant  14: 

9(3)  62474.3925(3+2.1700) 

«*(3)  “  (3  +  1.3559  ±i2.9567)(3+  10.8140  +  j7.7387)(s+  100) 

Plant  15: 

9(3)  202515.9298(3  +  4.7795) 

S.(s)  ~  (3  +  2.7986  ±;5.2958)(3+  10.8140  + j7.7387)(3+  100) 

Plant  16: 

9(3)  137680.4259(3  +  3.0148) 

«e(s)  “  (3  +  1.7875  +  i4.0287)(3  +  10.8140  +  j7.7387)(3  +  100) 
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C.t  Longitudinal  P4(s)  Short-Period  Transfer  Functions 


Plant  1: 


Plant  2: 


Plant  3: 


Plant  4: 


Plant  5: 


Plant  6: 


Plant  7: 


Plant  8: 


Plant  9: 


Plant  10: 


Plant  11: 


9(f) 

d(s) 


=  0 


9(f) 

dis) 


=  0 


9(^) 

dis) 


=  0 


9(f) 

dis) 


=  0 


qis)  2.0708(5  +  3.1857) 
rf(5)  “  (a  +  1.8419  ±j4.4824) 


^(5)  _  1.5778(5  +  2.1433) 

dis)  ~  (5  +  1.2778  +  ^3.9504) 


qjs)  _  5.1162(5  +  4.4250) 

dis)  '  (5  + 2.6377  ±j7.0929) 
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Plaat  12: 


Plant  13: 


Plant  14: 


Plant  15: 


Plant  16: 


q(a)  _  3.89€2(s  +  3.3681) 

d(s)  ~  (a  +  2.0079  ±  j6.2078) 

g(a)  _  2.0071(5-1-3.0415) 

d(5)  “  (a  -I- 1.7809  ±  j3.3700) 

q(a)  1.5290(s -I- 2.1700) 
dU)  ~  («  +  1-3559  ±  i2.9567) 

q(a)  _  4.9563(a-f  4.7795) 

rf(a)  “  (a  -»-  2.7986  ±  j5.2958) 
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Appendix  D.  Lateral  Plant  Transfer  Functions 


This  appendix  lists  the  effective  plant  transfer  functions  P(„(<s)  and  the  disturbance  transfer 
functions  Ps^tis)  for  the  lateral  channel. 


D.l  Transfer  Function  P«„(3)  = 


Plant  1: 

p(s)  160.79s(s  +  0.6633  ±  j2.5958)(s  +  4.4699)(s  +  1968.5371) 

S„(s)  ~  (s-  0.0226)(s  +  4.4372)(s  +  4.7342)(s  +  100)(a  +  0.6184  ±  i2.5351)(a  +  11.0737  ±  jlO.0893) 


Plant  2: 

p(s)  80.39s(s  +  0.5265  i:j2.3165)(s  +  4.4466)(s  + 7793.0854) 

S„(s)  ~  (s~  0.0193)(s  +  4.4372)(s  +  4.7718)(s  +  100)(s  +  0.5299  ±  j2.3482)(s  +  11.0737  ±  ilO.0893) 


Plant  3: 

p(s)  122.33s(s  +  0.4805  ±  j2.2153)(s  +  4.4576)(s  +  3908.3631) 

«,(a)  “  (s  ~  0.0219)(«  +  4.4372)(a  +  3.7251)(«  +  100)(«  +  0.4638  ±  j2.2221){s  +  11.0737  ±  jlO.0893) 


Plant  4: 

p{s)  122.33s(s  +  0.5062  ±  j2.2720){s  +  4.4777)(«  +  1966.3849) 

S,{s)  ~  {s-  0.0225)(s  +  4.4372)(s  +  3.6224)(s  +  100)(s  +  0.4629  ±  j2.2231)(5  +  11.0737  ±  jlO.0893) 


Plant  5: 

p{s)  61.17s(s  +  0.4009  ±j2.0272){s+4.4485)(s  + 7791.2275) 

««(s)  ~  (a  -  0.0193)(a  +  4.4372)(s  +  3.6473)(a  +  100)(5  +  0.3970  ±  j2.0586)(s  +  11.0737  ±  jlO.0893) 


Plant  6; 

p(s) _ 61.17a(a  +  0.4121  ±  j2.0544)(a  +  4.4596)(a  +  3907.3803) _ 

S„{s)  ~  (s-  0.0198)(s  +  4.4372)(a  +  3.5443)(s  +  100)(a  +  0.3962  ±  j2.0595)(5  +  11.0737  ±  jlO.0893) 
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Plant  7: 

jKf).  _ _  397.05a{s  +  0.9717  ±  >3.9639)(s  +  4.4840)(s  +  3915.6867) 

tf.(a)  ”  («  -  0.0142)(«  +  4.4372)(3  +  7.6264)(s  +  100)(«  +  0.9762  ±  j3.9567)(s  +  11.0737  ±  jlO.0893) 

Plant  8: 

p(a)  397.05a(«  +  1.0049  ±  j4.Q547)is  +  4.5286)(5  +  1973.6692) 

^.(«)  ~  -  0.0146)(a  +  4.4372)(a  +  7.4112)(«  +  100)(»  -|-  0.9755 ±  j3.9582)(«  +  11.0737  ±  jlO.0893) 


Plant  9: 

p(3)  _  198.52i(«  +  0.8191  ±  j3.6352)(a  +  4.4599)(s  +  7797.5219) 

^a(«)  ~  -  0.0125)(d  +  4.4372)(a  +  7A753){s  +  100)(5  +  0.8353  ±  j3.6669)(5  +  11.0737  ±  j 10.0893) 


Plant  10: 

p(a)  _  302.09s(s  +  0.7404  ±  j3A732)(a  +  4.4808)(a  +  3912.2895) 

f,(a)  ~  (a-  0.0142)(«  +  4.4372)(a  +  5.8217)(s  +  100)(5  +  0.7348  ±  j3A662)(a  +  11.0737  ±  jlO.0893) 

Plant  11: 

p(«)  151.04a(«  +  0.6236  ±  i3.1825)(«  +  4.4585)(3  +  7794.6027) 

««(«)  ~  -  0.0125)(a  +  4.4372)(3  +  5.7034 )(a  +  100)(s  +  0.6291  ±  j3.2110)(ii  +  11.0737  ±  j  10.0893) 

Plant  12: 

p{a)  _  265.76«(s  +  0.8352  ±  j3.3289)(«  +  4.4962)(5  +  1971.1031) 

«„(«)  ~  {a-  0.0177)(«  +  4.4372)(s  +  6.0711)(s  A  100)(«  +  0.7970  ±  j3.2454)(s  +  11.0737  ±  jlO.0893) 


Plant  13: 

f(3) _ 202.203(3  +  0.6111  ±  j2.8454)(a  +  4.4675)(3  -f  3910.3260) _ 

«a(3)  ”  (3  -  0.0173)(3  +  4.4372)(3  +  4.7717)(3  +  100)(3  +  0.5996  ±  j2.8429)(3  +  11.0737  ±  jlO.0893) 


D-2 


D,t  Tromftr  Function  P,„(3)  = 


Plant  1: 

p(«)  160.79«(j  +  5.1399  ±  >2.9560) 

Sris)  ~  («  -  0.0226)(3  +  4.4372)(a  +  4.7342)(«  +  100)(«  +  0.6184  ±  >2.5351) 

Plant  2: 

p(a) _ 80.39a(a  +  4.4056  ±  >3.2242) _ 

~  (s-  0.0193)(a  +  4.4372)(i  +  4.7718)(s  +  100)(«  +  0.5299  ±  >2.3482) 

Plant  3: 

p(a) _ 122.33a(a  +  3.9106  ±  >3.3115) _ 

Sris)  ~  0.0219)(a  +  4.4372)(«  +  3.7251)(«  +  100)(5  +  0.4638  ±  >2.2221) 

Plant  4: 

p(s)  122.33«(«  + 3.9106  ±>3.3115) 

«r(5)  ~  (a  -  0.0225)(a  +  4.4372)(a  +  3.6224)(«  +  100)(«  +  0.4629  ±  >2.2231) 

Plant  5: 

p(a)  61.17«(«  + 3.3519  ±>3.3289 

tfr(«)  ~  («  -  0.0193)(s  +  4.4372)(a  +  3.6473)(«  +  100)(«  +  0.3970  ±  >2.0586) 

Plant  6: 

p{s)  61.17«(«  + 3.3519  ±>3.3289) 

Sris)  ~  (a  -  0.0198)(a  +  4.4372)(a  +  3.5443)(a  +  100)(a  +  0.3962  ±  >2.0595) 

Plant  7: 


pis)  397.04a(a  + 8.0769  ±>4.6451) 

S^  ~  is-  0.0142)(a  +  4.4372)(a  +  7.6264)(a  +  100)(a  +  0.9762 ± >3.9567) 
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Plant  8: 


|>(a)  _ _ 397.05a(a  4  8.0769  ±  >4.6451) _ 

Sr{i)  ~  {s-  0.0146)(3  +  4.4372)(3  +  7.4112)(5  +  100)(3  +  0.9755  ±  >3.9582) 

Plant  9: 

pjs) _ 198.52a(a  4-  6.9231  ±  >5.0666) _ 

«,(a)  “  (a  -  0.0125)(a  +  4.4372)(a  +  7.4753)(a  +  100)(a  +  0.8353  ±  >3.6669) 

Plant  10: 

p(a)  302.09a(a  +  6.1452  ±  >5.2038) 

^a)  “  (a  -  0.0142)(a  +  4.4372)(a  +  5.8217)(a  +  100)(a  +  0.7348  ±  >3.4662) 

Plant  11: 

p(a) _ 151.05a(a  4-  5.2673  ±  >5.2311) _ 

Sr(s)  “  (a  -  0.0125)(a  +  4.4372)(a  +  5.7034)(a  +  100)(a  +  0.6291  ±  >3.2110) 

Plant  12: 

p(a)  265.76a(a  + 6.6080  ±>3.8003) 

«r(a)  “  (a  -  0.0177)(a  +  4.4372)(a  +  6.0711)(a  +  100)(a  +  0.7970  ±  >3.2454) 

Plant  13: 

p{s)  202.20a(a  + 5.0275  ±>4.2574) 

«,(a)  “  (a  -  0.0173)(a  +  4.4372)(a  +  4.7717)(a  +  100)(a  +  0.5996  ±  >2.8429) 
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D.3  Trnntfer  Function  P,„(«)  = 


Plant  1: 


$(s)  49.83(«  -  0.4498)(  J  -  33.9420)(a  +  0.7.9419  ±  j7.9450)(s  +  5.0220) 

Sr{»)  ~  -  0.0226)(3  +  4.4372)(a  +  4.7342)(a  +  100)(a  +  0.6184  i  j2.5351)(«  +  11.0737  ±  jlO.0893) 


Plant  2: 


${o)  24.92(a  +  0.7190)(«  +  12.6620)(a  +  15.5574  ±  j22.4002){a  +  4.0981) 

Sr{s)  ”  (a  -  0.0193)(a  +  4.4372)(a  +  4.7718)(a  +  100)(a  +  0.5299  ±  j2.M82){»  +  11.0737  ±  >10.0893) 


Plant  3: 


/3(a)  37.91(a  -  4.3592)(a  -  25.1150)(a  +  4.3592)(a  +  6.4502  ±  >5.1924) 

Sris)  ~  (a  -  0.0219)(a  -|-  4.4372)(a  +  3.7251)(a  +  100)(a  +  0.4638  ±  >2.2221)(a  +  11.0737  ±  >10.0893) 


Plant  4: 


^(a) _ 37.91(a  -  0.4711)(a  -  34.5206)(a  +  3.0440)(a  -j-  8.5081  ±  >8.7367) _ 

Sris)  ’  (a  -  0.0225)(a  +  4.4372)(a  +  3.6224)(a  +  100)(«  +  0.4629  ±  >2.2231)(a  +  11,0737  ±  >10.0893) 


Plant  5: 

/3(a)  18.96(a  +  0.5629)(a  +  6.0674)(a  +  8.9243)(a  -  15.3000  ±  >21.4138) 

«r(a)  ~  (a  -  0.0193)(a  +  4.4372)(a  +  3.6473)(a  +  100)(a  +  0.3970  ±  >2.0586)(a  +  11.0737  ±  >10.0893) 

Plant  6: 

/3(a)  18.96(a  -  3.6589)(a  -  25.7908)(a  +  1.5782)(a  +  6.3605  ±  >5.9201) 

^r(3)  “  («  -  0.0198)(a  +  4.4372)(a  +  3.5443)(a  +  100)(a  +  0.3962  ±  >2.0595)(a  +  11.0737  ±  >10.0893) 


Plant  7: 

/3(a)  78.31(a  -  1.5107)(a  -  48.1421)(a  +  15.3632)(a  +  0.3294  ±  >4.9133) 

tf,(a)  "  (a  -  0.0142)(a  +  4.4372)(a  +  7.6264)(a  +  100)(a  +  0.9762  ±  >3.9567)(a  +  11.0737  ±  >10.0893) 
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Plant  8: 


0{s)  78.31(a  -  0.2450)(a  -  56.7801)(a  +  11.8098)(«  +  5.6836  ±  j7.7627) 

Sria)  ~  (a-0.0146)(«  +  4.4372)(«  +  7.4112)(«+  100)(«  + 0.9755  ±j3.9582)(5 -Ml. 0737 ±>10.0893) 

Plant  9: 

fi(s)  39.15(3  =  0.9761)(3  +  2.2689)(3  ±  17.0154){3  -  27.1719  ±  >20.4655) 

Sria)  ~  {a-  0.0125)(3  +  4.4372)(3  +  7.4753)(3  +  100)(3  +  0.8353  ±  >3.6669)(3  ±  11.0737  ±  >10.0893) 

Plant  10: 

0ia)  59.58(3  -  1.5089)(3  -  49.4974)(3  +  11.2967)(3  +  1.8778  ±  >4.5724) 

~  (3  -  0.0142)(3  +  4.4372)(3  +  5.8217)(3  +  100)(3  +  0.7348  ±  >3.4662)(3  +  11.0737  ±  >10.0893) 

11: 

29.79(3  +  0.7439)(3  +  3.1626)(3  +  13.4091)(3  -  26.8070  ±  >18.2492) 

~  (3  -  0.0125)(3  +  4.4372)(3  +  5.7034)(3  +  100)(s  +  0.6291  ±  >3.2110)(3  +  11.0737  ±  >10.0893) 

Plant  12: 

0(3)  64.07(3  -  0.3155)(3  -  45.3608)(3  +  8.4734)(3  +  6.7694  ±  >7.5346) 

^r(3)  ~  (9-  0.0177)(3  +  4.4372)(3  +  6.0711)(3  +  100)(3  +  0.7970  ±  >3.2454)(3  +  11.0737  ±  >10.0893) 

Plant  13: 

0ia)  48.74(3  -  2.4083)(3  -  37.6040)(3  +  7.5162)(3  +  3.5545  ±>3.4846) 

«r(a)  ~  {a-  0.0173)(3  ±  4.4372)(3  +  4.7717)(3  +  100)(3  ±  0.5996  ±  >2.8429)(3  +  11.0737  ±  >10.0893) 


Sria) 

Plant 

^is) 
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D.4  Transfer  Function  Pe„(3)  = 


Plant  1: 

/3(«)  49.83(«  -  0.0494)(i  -  40.3037)(a  +  4.7197) 

tfr(a)  ~  {s-  0.0226)(«  +  4.4372)(«  +  4.7342)(3  +  100)(a  +  0.6184  ±  j2.5351) 

Plant  2: 

0(s)  24.92(3  -  0.0419)(3  -  40.4971)(3  +  4.7557) 

Sr{s)  ~  (s-  0.0193)(3  +  4.4372)(3  +  4.7718)(3  +  100)(3  +  0.5299  ±  >2.3482) 

Plant  3: 

0is)  37.91(3  -  0.0476)(3  -  40.6275)(3  +  3.7015) 

Sr(s)  ~  (3  -  0.0219)(3  +  4.4372)(3  +  3.7251)(3  +  100)(3  +  0.4638  ±  >2.2221) 

Plant  4: 

^(s)  37.9143(3  -  0.0490)(3  -  40.6275)(3  +  3.5978) 

Sr(s)  ~  (s-  0.0225)(3  +  4.4372)(3  +  3.6224)(3  +  100)(3  +  0.4629  ±  >2.2231) 

Plant  5: 

/J(s)  _  18.96(3 -0.0415)(s-40.7744)(3  +  3.6233) 

^r(3)  ~  (s-  0.0193)(3  +  4.4372)(3  +  3.6473)(s  +  100)(3  +  0.3970  ±  >2.0586) 

Plant  6: 

/3(s)  18.96(3  -  0.0427)(3  -  40.7745)(3  +  3.5194) 

«/(3)  ~  (3  -  0.0198)(3  +  4.4372)(3  +  3.5443)(3  +  100)(3  +  0.3962  ±  >2.0595) 

Plant  7: 

j3(s)  78.31(3  -  0.0307)(3  -  63.3574)(3  +  7.6099) 

«r(3)  ”  (3  -  0.0142)(3  +  4.4372)(s  +  7.6264)(3  +  100)(3  +  0.9762  ±  >3.9567) 
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Plsuit  8: 

$(s)  78.31(s-0.0316)(5-63.3575)(5+  7.3937) 

^r(s)  ~  («  -  0.0146)(s  +  4.4372)(5  +  7.4112)(3  +  100)(5  +  0.9755  ±  j3.9582) 

Plant  9: 

I3{s)  39.15(a  -  63.6605)(s  -  0.0267)(s  +  7.4565) 

rfr(3)  ~  («  -  0.0125)(s  +  4.4372)(«  +  7.4753 )(«  +  100)(«  +  0.8353  ±  >3.6669) 

Plant  10: 

/3(s)  59.58(s  -  0.0304)(5  -  63.8655)(3  +  5.7945) 

«,(3)  ~  (a  -  0.0142)(a  +  4.4372)(s  +  5.8217)(a  +  100)(«  +  0.7348  ±  >3.4662) 

Plant  11: 

/3(s) _ 29.79(3  -  0.0265)(8  -  64.0958)(a  +  5.6767) _ 

rf,(3)  ~  (s-  0.0125)(3  +  4.4372)(3  +  5.7034)(3  +  100)(3  +  0.6291  ±  >3.2110) 

Plant  12: 

)9(3)  64.07(3  -  0.0385)(3  -  51.8282)(3  +  6.0553) 

^r(«)  ”  (a  -  0.0177)(3  +  4.4372)(3  +  6.0711)(3  +  100)(3  +  0.7970  ±  >3.2454) 

Plaint  13: 

/3(s)  48.74(3  -  0.0371)(3  +  52.2440)(3  +  4.7468) 

^r(3)  “  (a  -  0.0173)(3  +  4.4372)(s  +  4.7717)(s  +  100)(j>  +  0.5996  ±  >2.8429) 
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D.S  Tmnsfer  Function  Prf„(s)  =  ^ 


Plant  1; 


Plant  2; 


Plant  3; 


Plant  4: 


Plant  5: 


Plant  6: 


Plant  7: 


Plant  8: 


Plant  9: 


Plant  10: 


Plant  11: 


p(«)  _  2.40a(a  +  0.6182  ±  i2.4803) 

d(s)  ~  (a  -  0.0226)(a  +  4.7342)(3  +  0.6184  ±  j2.5351) 

p(s)  2.40a(a  +  0.5299  ±j2.3023) 

d{a)  ~  (s-  0.0193)(s  +  4.7718)(a  +  0.5299  ±  j2.3482) 

p{a)  1.83«(a  +  0.4703  ±j2.1729) 

~  (a  -  0.0219)(a  +  3.7251)(3  +  0.4638  ±  j2.2221) 

p{a) _ 1.83a(a  +  0.4703  ±  j2.1729 _ 

d(a)  ~  (s-  0.0225)(a  +  3.6224)(a  +  0.4629  ±  j2.2231) 

p(3)  1.83a(a  +  0.4031  ±j2.0156) 

d(a)  ~  (a  -  0.0193)(a  +  3.6473)(a  +  0.3970  ±  ;2.0586) 

p(a)  1.83a(«  +  0.4031  ±  j2.0156) 

d(a)  ”  (a  -  0.0198)(a  +  3.5443)(a  +  0.3962  ±  j2.0595) 

p{3)  5.93a(a  +  0.9714  ±  i3.8977) 

d(a)  ~  (a  -  0.0142)(a  +  7.6264)(a  +  0.9762  ±  j3.9567) 

j>(a)  _  5.93a(a  +  0.9714  ±  j3.8977) 

d(a)  ~  (a  -  0.0146)(a  +  7.4112)(a  +  0.9755  ±  ;3.9582) 

p(a)  5.93a(a  +  0.8326  ±  j2.6179) 

d(a)  ~  (a  -  0.0125)(a  +  7.4753)(a  +  0.8353  ±  >3.6669) 

p(a) _ 4.51a(a  +  0.7391  ±  >3.4145) _ 

d(a)  ~  (a  -  0.0142)(a  +  5.8217)(a  +  0.7348  ±  >3.4662) 

p(a)  _ _ 4.51a(a  +  0.6335  ±  >3.1674) _ 

d(a)  ~  (a  -  0.0125)(a  +  5.7034)(a  +  0.6291  ±  >3.2110) 
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Plant  12: 


p{s)  _  _ 3.973(3  +  0.7947  ±  J3.1888) _ 

d(s)  ~  (3  -  0.0177)(3  +  6.0711)(s  +  0.7970  ±  i3.2454) 

Plant  13: 

pjs) _ 3.023(3  +  0.6047  ±  j2.7935) _ 

d{s)  ~  (3-0.0173)(3  +  4.7717)(3  4  0.5996  ±  j2.8429) 
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D.6  Transfer  Function  Pd„(a)  =  ^ 


Plant  1: 


Plant  2: 


Plant  3: 


Plant  4: 


Plant  5: 


Plant  6: 


Plant  7: 


Plant  8: 


Plant  9: 


Plant  10: 


Plant  11: 


0{s)  0.0346(s  +  0.5361)(a  -  36.6976) 

dls)  ~  (s-  0.0226)(s  +  4.7342)(5  +  0.6184  ±  j2.5351) 

/?(s)  0.0346(a  -  1.9495)(«  -  34.1827) 

d(s)  ~  {a-  0.0193)(3  +  4.7718)(«  +  0.5299  ±  j2.3482) 

0{s)  0.0264(3  +  0.4604)(3  -  32.4824) 

d(s)  ~  (3 -0.0219)(3  +  3.7251)(s+ 0.4638  ±j2.2221) 

^{a)  _  0.0264(3  +  0.4604)(3  -  32.4824) 

d(a)  ~  (3 -0.0225)(3  +  3.6224)(3  + 0.4629  ±j2.2231) 

d{a)  0.0264(3  +  0.4193)(3  -  30.5580) 

d{a)  ~  (3 -0.0193)(3  +  3.6473)(3  + 0.3970  ±j2.0586) 

/3(3)  0.0264(3  +  0.4193)(3  -  30.5578) 

d(a)  “  (3  -  0.0198)(3  +  3.5443)(3  +  0.3962  ±  i2.0595) 

/9(a)  0.0855(3  +  0.4961)(3  -  39.6594) 

d(a)  ~  (3 -0.0142)(3  +  7.6264)(3  + 0.9762  ±j3.9567) 

^(3)  _  0.0855(3  +  0.4961)(3  -  39.6595) 

d(a)  ~  (a-  0.0146)(3  +  7.4112)(s  +  0.9755 ± j3.9582) 

J3(a)  0.0855(3  +  0.4713)(s  -  35.7606) 

d(s)  ~  (3  -  0.0125)(3  +  7.4753)(3  +  0.8353  ±  j3.6669) 

/3(a)  0.0855(3  +  0.4514)(s  -  33.1269) 

d(a)  ~  (a-  0.0142)(s  +  5.8217)(s  +  0.7348  ±  j3.4662) 

/3(a)  0.0651(3 +  0.4250)(3- 30.1487) 

d(a)  ~  (a-  0.0125)(s  +  5.7034)(3  +  0.6291  ±  j3.2110) 
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Plant  12: 


P(a) _ 0.0572(3  -f  0.5257)(^  -  37.4245) 

d{s)  ~  (3 -0.0177)(3  +  6.0711)(3  + 0.7970  ±  >3.2454) 

Plant  13: 

^{s) _ 0.0436(3  +  0.4666)(3  ~  32.0527) 

d(3)  ”  (3  -  0.0173)(3  +  4.7717)(3  +  0.5996  ±  >2.8429) 
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Appendix  E.  MISO  Equivalent  Transfer  Functions 

This  appe&dix  lists  the  MISO  equivalent  transfer  functions  used  for  roll  and  yaw  channel 
loop-shaping. 

E.l  Roll  Channel  q2i{s) 

Plant  1: 

.  .  _ _ 312234.51s(s  -  40.5188) _ 

^  (s  -  0.0494)(s  -  40.3037)(s  4.7197)(s  +  100)(s  -f- 11.07365  ±  j  10.0893) 

Plant  2: 

^  ~  («  -  0.0494)(a  -  40.4971)(j  -|-  4.7557)(5  -|-  100)(j(  -|- 11.07365  ±  jlO.0893) 

Plant  3: 

/X  _ _ 475115.98s(s  -  40.8274) _ 

”  (a  -  6.0476)(5  -  40.6275)(3-f-  3.7015)(5  -|-  100)(s  +  11.07365  ±  >10.0893) 

Plant  4: 

.X  _ _ 237557.98s(s  -  40.8274) _ 

*  (s  -  0.0494)(5  -  40.6275)(s  +  3.5978)(s  -{■  100)(s  -|- 11.07365  ±  jlO.0893) 

Plant  5: 

^  («  -  0.0415)(s  -  40.7744)(s  +  3.6233)(s  -|-  100)(s  -|- 11.07365  ±  jlO.0893) 

Plant  6; 

/X  _ _ 237557.983(s  -  40.9676) _ 

^  («  -  0.0427)(s  -  40.7745)(s  +  3.5194)(s  +  100)(s  -|- 11.07365  ±  ;10.0893) 
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Plant  7: 


.  V  _ _ 1542055.60^(3  -  63.6723) _ 

ftais;  -  _  o.o307)(3  -  63.3574)(3  +  7.6099)(3  +  100)(«  +  11.07365  ±  j  10.0893) 

Plant  8: 

.  . _ 771027.793(3  -  63.6723) _ 

q72[s)  -  o.o316)(3  -  63.3575)(3  +  7.3937)(3  +  100)(3  +  11.07365  ±  jlO.0893) 

Plant  9; 

,  .  _  _ 1542055.603(3  -  63.9620) _ 

-  (3  _  o.0267)(3  -  63.6605)(3  +  7.4565)(3  +  100)(3  +  11.07365  ±  >10.0893) 

Plant  10: 

^  _  1173245.303(3  -  64.1573) 

(a  _  0.0304)(3  -  63.8655)(3  +  5.7945)(3  +  100)(3  +  II .  ; 7365  ±  >10.0893) 

Plant  11; 

,  .  _ _ 1173245.303(3  -  64.3777) _ 

(3  -  0.0265)(3  -  64.0958)(3  +  5.6767)(3  +  100)(3  +  11.07365  ±  >10.0893) 

Plant  12: 

,  V  _  517074.413(3  -  52.0921) 

(3  -  0.0385)(3  -  51.8282)(3  +  6.0553)(3  +  100)(3  +  11.07365  ±  >10.0893) 

Plant  13: 

,  .  _ _ 785291.763(3  -  52.4889) _ 

^  ~  (3  -  0.0371)(3  -  52.2440)(3  +  4.7468)(3  +  100)(3  +  11.07365  ±  >10.0893) 
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Plant  1: 


63417.73(«  -  40.3100)(s  -  40.5188) 

****^^'  ~  (s  -  40.2054)(s  +  4.4706 )(a  +  1293.9194)(a  +  100)(a  +  0.6847  ±  j2.5895) 

Plant  2: 

132266.67(«  -  40.5086)(tf  -  40.7031) 

**®*^^^  “  (s  -  40.4193)(a  +  4.4466)(a  +  5342.83)(«  +  100)(a  +  0.5219  ±  i2.3163) 

Plant  3: 

97956.56(a  -  40.6363)(«  -  40.8274) 

"  (a  -  40.5431)(«  +  4.4553)(«  +  2613.6590)(a  +  100)(s  +  0.4794  ±  j2.2127) 

Plant  4: 

. _ 48912.58(3  -  40.6321)(fl  -  40.8274) _ 

?33«(«)  -  _  4o.5393)(5  ^  4.4629)(3  +  1318.9502Xs  +  100)(3  +  0.5237  ±  j2.2633) 

Plant  5; 

97943.79(3  -  40,7829)(3  -  40.9676) 

~  (3  -  40.6915)(3  +  4.4470)(3  +  5201.9369)(3  +  100)(3  +  0.3964  ±  j2.0268) 

Plant  6: 

.  _  48799(3  -  40.7788)(3  -  40.9676) 

**®*^*'  “  (3  -  40.6867)(3  +  4.4501 )(3  +  2612.6786)(3  +  100)(3  +  0.4117  ±  j2.0498) 

Plant  7: 

203390.43(3  -  63.3736)(3  -  63.6723) 

933e(«)  -  _  63.2294)(3  +  4.4885)(3  +  2620.9281)(3+  100)(3  +  0.9724  ±  j3.9645) 
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Plant  8: 


9ss«(«)  =  77 

Plant  9: 

9a3*{s)  =  ^ 

Plant  10: 

154544.77(a  -  63.877l)(s  -  64.1573) 

~  (s~  63.7418)(j  +  4.4832)(tf  +  2617.5386)(5  +  100)(s  +  0.7413  ±  j3.4726) 

Plant  11: 

.  .  _  154212.1633(«  -  64.1069)(s  -  64.3777) 

**®*^^^  “  (a  -  63.9738)(3  +  4.4597)(3  +  5205.2705)(s  +  100)(jj  +  0.6204  ±  j3.1827) 

Plant  12: 

.  .  83727.5912(«  -  51.8356)(s  -  52.0921) 

~  (s  -  51.7151)(s  +  4.5037)(«  +  1323.6283)(3+  100)(«  +  0.8576  ±  j3.3283) 

Plant  13: 

.  ,  _ 126250.73(a  -  52.2544)(5  -  52.4889) _ 

~  (5  -  52.1406)(s  +  4.4682)(5  +  2615.5985)(3  +  100)(s  +  0.6115  ±  j2.8439) 


102593.59(5  -  63.3658)(5  -  63.6723) 

-  63.2232)(3  +  4.5416)(5+  1326.1636)(5+  100)(5+  1.0261  ±  >4.0568) 


_ 202804.46(5  -  63.6758)(5  -  63.9620) _ 

-  63.5351)(5  +  4.4620)(5  +  S208.1863)(5  +  100)(5  +  0.8152  ±  >3.6363) 


E-4 


Appendix  F.  Difference  Equations 


This  app«idix  presents  the  difference  equations  needed  for  digital  implementation  of  the 
designed  flight  controllers  and  prefilters.  The  equations  assume  a  sampling  rate  of  50Hz. 

Equations  (4.15)  and  (4.13)  are  the  s-plane  representations  of  the  designed  longitudinal 
prefilter  and  controller,  respectively.  Equations  (5.19)  and  (5.20)  are  2-plane  representations 
of  the  lateral  prefilters.  Equations  (5.15)  and  (5.16)  are  z-plane  representations  of  the  lateral 
contrdlers.  The  prefilters  and  controllers  are  converted  to  difference  equations  for  implementation 
on  Lambda’s  digital  flight  controller. 

The  z-plane  transfer  functions  are  of  the  form: 

Uiz)  d„z"-J-d„-i-?’‘-‘-l-...  +  d)^  +  d,  ^ 

where:  n  >  m  and  dn  =  !•  The  difference  equation  is  formed  by  dividing  the  equation  by  z",  cross 
multiplying,  and  replacing  the  negative  powers  of  z  with  the  appropriate  delay  to  obtain: 

vikT)  =  -d„.Mik-l)T]-...-dty[ik  +  l-n)T]-doylik-n)T] 

-|-c,„u[(fc  -I-  m  -  n)r]  -I-  c,„_iu[(*r  +  m-n  -  1)T]  +  ... 
-|-Citt[(fc-|-l-n)r]-|-Co«[(*:-n)T]  (F.2) 

The  difference  equations  are  listed  below: 

fu(kT)  =  2,57754514/i,[(ib  -  l)r] 

-2.20234613/„[(jfc  -  2)T] 

-0.62346973/ii[(ifc  -  3)r] 

-|-0.03569841gre/[ikT] 

-0.03436716g„,[(Jt  -  1)T]  (F.3) 

Pu(fcT)  =  1.46358316pu[(fc  -  1)T] 

-0.028495385ii[(jfe  -  2)r] 
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-0.34868264^„{(Jb  -  3)7] 

-0.08121895^1, ((ib  -  4)7] 

-0.005186185u[(jb  -  5)7] 

+12.25155696^,,.,[jfc7] 

-29.6970390(We«,[(ib-  1)7] 

+11.84532698rfe..J(ifc  -  2)7] 

+23.14265795^.„,[(ib  -  3)7] 

-24.09492395^.„,[(ik  -  4)7] 

+6.55634105^...,[(ib  -  5)7]  (F.4) 


fn{kT)  =  1.70497427/„[(jk  -  1)7] 

-0.72555746/„[(ib  -  2)7] 

0.02058319p,./[fc7]  (F.5) 


gn(kT)  =  1. 5584541 lff„[(Jk  -  1)7] 

-0.17971014^aj{(fc  -  2)7] 

-0.318840589a,[(jfe  -  3)7] 

+0.05700483^ja[(Jfc  -  4)7] 

+0.00289855ffjal(ib  -  5)7] 

+5.99046148^..^[ifc7] 

-15.191260285.,,,[(ib-  1)7] 

+7.20757487«..^,[(jk  -  2)7] 

11.23289141«.„,[(jk  -  3)7] 

- 13.19796181#..^  J(jk  -  4)7] 

+3.95844341#.„.[(fc  -  5)7]  (F.6) 


fsaikT)  =  2.90477497/33[(Jfc-l)7] 

-2.81330221/33[(jk  -  2)7] 

-0,90846623/33[(/fc  -  3)7] 

+0.00308112/?re/[ik7] 

-0.00302011/?„;[(ifc  -  1)7]  (F.7) 
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i733(*r) 


3.44977611s33{(it  -  1)T] 
-4.551822995ss[(ib  -  2)T] 
+2.84736852533((*  -  3)T] 
-0.838S8820533l(*  -  4)7’] 
+0.09326655ff3s[(fc  -  5)r] 
-22.12114283^r..,l*T] 
+63.99586422«r«J(ifc  -  1)T] 
-39.63341388«r..4l(ik  -  2)7^ 
-44.11914690^,..,[(ib  -  3)r] 
+61.75453156rfr..,[(*  -  4)7’] 
-19.87674225^r.«,[(Jt  -  5)7’] 


(F.8) 
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